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ABSTRACT 

We report the results of the counterpart identification and a detailed analysis of the 
physical properties of the 48 sources discovered in our deep 1.1mm wavelength imag- 
ing survey of the GOODS-South field using the AzTEC instrument on the Atacama 
Submillimeter Telescope Experiment (ASTE). One or more robust or tentative coun- 
terpart candidate is found for 27 and 14 AzTEC sources, respectively, by employing 
deep radio continuum, Spitzer MIPS & IRAC, and LABOCA 870 iim data. Five of 
the sources (10%) have two robust counterparts each, supporting the idea that these 
galaxies are strongly clustered and/or heavily confused. Photometric redshifts and 
star formation rates (SFRs) are derived by analyzing UV-to-optical and IR-to-radio 
SEDs. The median redshift of z me d ~ 2.6 is similar to other earlier estimates, but we 
show that 80% of the AzTEC-GOODS sources are at z ^ 2, with a significant high 
redshift tail (20% at z > 3.3). Rest-frame UV and optical properties of AzTEC sources 
are extremely diverse, spanning 10 magnitude in the i— and A"— band photometry (a 
factor of 10 4 in flux density) with median values of i — 25.3 and K = 22.6 and a 
broad range of red colour (i — K =0-6) with an average value of i — K k, 3. These 
AzTEC sources are some of the most luminous galaxies in the rest-frame optical bands 
at z ^ 2, with inferred stellar masses = (1-30) xlO lo M0 and UV-derived star for- 
mation rates of SFR UV > IO^Mq yr" 1 . The IR-derived SFR, 200-2000 M yr" 1 , 
is independent of z or A/*. The resulting specific star formation rates, SSFR w 1-100 
Gyr -1 , are 10-100 times higher than similar mass galaxies at z = 0, and they extend 
the previously observed rapid rise in the SSFR with redshift to z = 2 — 5. These 
galaxies have a SFR high enough to have built up their entire stellar mass within 
their Hubble time. We find only marginal evidence for an AGN contribution to the 
near-IR and mid-IR SEDs, even among the X-ray detected sources, and the derived 
M* and SFR show little dependence on the presence of an X-rav bright AGN. 
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1 INTRODUCTION 

Early studies of the far-infrared (FIR) cosmic background 
indicated that up to half of the cosmic energy density 
is ge nerated by dusty starbursts and active galactic nu- 
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revealed that this IR background is resolved into a large pop- 
ulation of discrete individual sources. 

Identifying and understanding the nature of these 
discrete FIR sources ( "submillimetre galaxies" or SMGs) 
has proven to be challenging because of the low angular 
resolution of these instruments and the faintness of counter- 
parts in the rest- frame optical and UV-bands (see review by 
iBIain et all 120021 ). Utilizing sub-arcsec astrometry of inter- 
ferometric radio continuu m data and sensit i ve spe ctroscopy 
using the Keck telescopes, IChapman et al.l (|2005h reported 
spectroscopic redshifts of 73 SMGs culled from earlier 
SCUBA surveys and concluded that they are massive, 
young objects seen during their formation epoch, with very 
high star formation rates at z > 1. Deep 24 [im band imag- 
ing usi ng the Mult i band Imaging Photometer for Spitzer 
(MIPS; iRieke et all |2004| ) on the Spitzer Space Telescope 
and f ollow-up spectrosco py using the Infrared Spectrograph 
(IRS: iHouck et aill2004h have also provided use ful insights 
on th e nature and redshifts of addi ti onal S M Gs llLutz et al 



Menendez-Delmestre et al.l 2007 ; Valiante et al 



Pope et all )2008al; iMenendez-Delmestre et all 120091: 



Huang et al] 120091 ). However, the use of high resolution 
radio continuum and MIPS 24 [im images for the coun- 
terpart identification suffers from a well-known systematic 
bias against high redshift (z > 3) sources. Indeed, a large 
fraction of the counterpart sources identified using direct 
interferometric imaging in the mm/submm wavelengths 
are shown to be extremely faint in nearly all other wave- 
length bands (r > 26, K > 24) with li t tle or no radio or 
MIPS 24 \im emission (llono et all 120061 ; IWang et all 120071 ; 
lYounger et ail 120071 . 120091 ), and high redshift SMGs may 
have been missed or mis-identified with a foreground source 
in the earlier studies. 

Obtaining a more complete understanding of the SMG 
population requires a study of a larger, more uniform sam- 
ple identified utilizing the deepest available multiwavelength 
complementary data and a robust counterpart identification 
method that is less prone to a redshift bias. In this pa- 
per we present the identification of the 48 AzTEC 1.1mm 
sources found in the deepest survey at mm wavelengths 
ever carried out in the Great Obs ervatories Origins Deep 
Survey-South (GOODS-S) field by IScott et all (|2010l . "Pa- 
per I" hereafter) . Several different identification methods are 
employed simultaneously to complement and to calibrate 
each other. A thorough analysis of the counterpart prop- 
erties and redshift distribution is also carried out as the 



GOODS-S field represents one of the most widely studied 
regions of sky with some of the deepest multi-wavelength 
data: X-ray data from Ch andra (|Luo et al.ll2008l ; IXue et all 
2011; Johnson et al.l [201 lh ■ optical to nea r-IR photometry 
from the Hubble Space Telescope (HST; iGiavalisco et alj 
l2004h . Spitzer IRAC (Dickinson et al. in prep.) and MIPS 
(Chary et al. in prep.) imaging in the mid-IR, submm 
imaging at 250 — 500 [im with the Balloon -borne Large 
Aper ture Submillimeter Telescope (BLAST; iDevlin et al.l 
120091 ) . and 1.4 GHz i nterferometric imaging with the Very 
Large Array (VLA; iKellermann et all 12003 ; iMiller et all 
l2008h . Extensive spectroscopy of optical sources in this field 
is also available dLe Fevre et al.l 20041; ISzokolv et al.l 2004 
Mignoli et"aill2005l;IVanzella et al.ll2005l . 120061: iNorris et al 
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2009; Trci ster et al.ll2009l; IWuvts et al.ll 2009; Bal estra et al 
201ol ; ISilverman et all |2010| ; ICasev et all l201ll). Including 
the AzTEC GOOD S-North field sources (|Perera et aill2008l ; 
IChapin et aT1l2009t ). our combined AzTEC-GOODS sample 
includes ~ 80 SMGs identified using a uniform set of cri- 
teria and the deepest multiwavelength data available and 
offers the best opportunity yet to examine the nature of 
the SMGs as a population and to verify the conclusions 
of earlier studies of most l y smaller and often radio-selected 
samples llLillv et al.|[l999l;lFox et al.ll2002l;lhdson et al.l | 2002| ; 



Webb et al.ll2003bl; IBorvs et al.ll2004l; IChapman et al.ll2005 
Ivison et al.ll2007l ; ]ciements et al.ll2008h - 



2 COUNTERPART IDENTIFICATION 

Here we describe the methods of identifying multiwavelength 
counterparts to the 48 Az TEC GOODS-S (A zTEC/GS here- 
after) sources reported bv lScott et al.l (|20ld ). We adopt the 
updated AzTEC source positions and photo metry derived 
using the improved point source kernel by iDownes et alj 
|201ll ). Our counterpart identification relies primarily on 
three observed multi- wavelength properties: (1) high reso- 
lution radio continuum; (2) Spitzer MIPS 24 um photome- 
try; and (3) red colors in the Spitzer IRAC bands. A robust 
counterpart is identified using a combination of these cri- 
teria for most AzTEC sources, and proposed identifications 
and multi-wavelength photometry for each of the AzTEC 
1.1mm sources are summarized in Table [T] and Tabled A 
more detailed discussion of the individual identification and 
the nature of the individual counterpart candidates are dis- 
cussed in Appendix [X] 



2.1 Methods 

Since the origin of the millimetre continuum emission de- 
tected by the AzTEC instrument is likely reprocessed ra- 
diation from dust-obscured starburst or AGN activity, the 
main data sets we examine for the multi- wavelength counter- 
part identification are mid-IR data from the Spitzer MIPS 
24 (j.m (full width at half maximum angular resolution of 
6>fwhm ~ 6") and IRAC 3.6 to 8.0 um band (6»fwhm ~ 2") 
and deep radio continuum data obtained using the VLA 
(#fwhm ~ 2"), exploiting the well-known r adio-IR correla - 
tion for star-forming galaxies (see review bv lCondon|[l99l ). 
The Spitzer IRAC and MIPS images and catalogues used 
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come from the Spitzer GOODIE the FIDElfl and the SIM- 
PLE0 Legacy Surveys. The radio continuum data used come 
from the VLA 1 . 4 GHz deep imaging su rvey {a ~ 8/^Jy; 
iMiller et ai1l2008l ; IKellermann et al.ll2008T ). Given their high 
resolution, the astrometric accuracy of these catalogues are 
sufficient to identify unique optical and near-IR counterparts 
in the deep ground-based telescope or HST images when 
such a counterpart is present. The i— band and K— band 
photometry of the counterpart candidates are also reported 
from the band-merged GOODS team photometry catalogue 
(Grogin et al., in prep.) cons tructed using a tem plate fit- 
ting software pa ckage TFIT (|Laidler et al.l 120071 ') and the 

alj~ 



MUSYC survey (|Gawiser et alj|2006h 



Unlike most previous works, we employed a variable 
search radius based on the beam size (Afwhm ~ 30") and 
the S/N of the AzTEC detection. Given the modest S/N 
(typically 5j 10), the positional offset between an AzTEC 
source and its counterpart is expected to be dominated by 
the map noise. This means we can exploit the measured S/N 
of each detection to constrain the counterpart identification. 
We derive the search radius, Rs, listed in Table [1] empiri- 
cally through simulations by injecting artificial sources into 
the signal map one at a time and measuring the distribution 
in the input to output source positions as described in Pa- 
per I. For each AzTEC source we select Rs such that there 
is a 95% probability that the true position of the source (as- 
sumed to be the position of the radio and/or Spitzer coun- 
terpart) will be within R s of the AzTEC centroid. 

The primary method of identifying AzTEC coun- 
terpar ts is the "P— statistic" described by iDownes et al.l 
(1986). This method computes the likelihood of a chance 
coincidence for each source in the comparison catalogue 
from the measured catalogue source density and the dis- 
tance to a given AzTEC source position. Following previ- 
ous work, a counterpart with a P— statistic less than 0.05 is 
deemed a "robust" identification, while a counterpart with 
0.05 < P < 0.20 is considered a "tentative" identification. 
Unlike most works, however, we compute all P-statistics 
based on the number density of all sources in the com- 
parison catalog, rather than the number density of sources 
brighter than the candidate counterpart in question. This 
means that all candidate counterparts equidistant from the 
AzTEC centroid will have the same P-statistic. This avoids 
biasing the identifications to the brightest radio and mid-IR 
sources, which could result in more misidentifications with 
low-redshift galaxies. 

For the radio data, we created a ^ 4a catalog using the 
SAD program in the Astronomical Image Processing Sys- 
tem (AlPSfl This program builds a catalog iteratively by 
first identifying bright pixels and then quantitatively test- 
ing their significance by fitting the PSF to the surrounding 
pixel brightness distribution. By allowing for a collection of 
connected sources as an acceptable model, this algorithm 
also provides a good estimate of the radio flux for extended 
objects as well. Submm/mm galaxies are almost always as- 
sociated with IRAC galaxies with faint but visible radio 



1 http: / /www. stsci.edu/science/goods/ 

2 http:/ /ssc. spitzer. caltcch.edu/legacy/abs/dickinson2. html 

3 http:/ /www. astro. yale.edu/dokkum/simple/ 

4 http://www.aips.nrao.edu/ 



emission. The extremely deep Spitzer data in these GOODS 
fields ensures that radio sour ces without an IRAC co unter- 
part are rare as reported by IKellermann et all (|2008l ). who 
find that only three out of 266 cataloged radio sources have 
no apparent counterpart at any other wavelengths. Taking 
advantage of this fact, we probe deeper into the radio data 
by creating a combined IRAC+VLA radio catalog by using 
the positions of IRAC sources detected with (S/N)s.e ^ 4 as 
prior positions. For each IRAC source, we fit a 2-D Gaussian 
to the radio map at the IRAC position, fixing the FWHM to 
1.6"and 2.8"in RA and Dec, respectivel y, based on the best - 
fit Gaussian to the synthesized beam (|Miller et al.l 12008). 
We limit the location of the peak to within 2"of the initial 
IRAC position. If the best-fit 1.4 GHz peak emission is > 3a 
of the rms noise in the surrounding region, we include this in 
our combined IRAC+VLA catalog. This list is cross-checked 
with the > 4a radio catalog created by the SAD program, 
and we use the SAD catalogue flux estimates where avail- 
able. The number density of IRAC+VLA sources in this 
catalog is 8330 deg" 2 for the GOODS-S+VLA catalog, and 
7860 deg~ 2 for the shallower SIMPLE+VLA catalog. 

For the MIPS 24 um catalogues, we use the number 
density of (S/N)24 1? 4 sources to compute the P-statistics, 
which are 45700 deg -2 and 25600 deg" 2 for the GOODS-S 
and FIDEL 24 urn catalogs, respectively. 

The third and an entirely new method we use for iden- 
tifying SMG counterparts employs their characteristic red 



millimetre continuum dlono et al. 


20061: IWang et all 12007. 


20 111; lYounger et al.l 2007. 200S 


, 200S 


• 


Hatsukadc et al. 


2010: Tamura et al. 2010: Ikarashi et al. 


1201 if) have shown 



that every source is detected in the IRAC 3.6 u.m and 4.5 urn 
bands at the ^ 1/iJy level, while their radio and MIPS 24 urn 
counterparts are not always detected in the best available 
data. By examining the spec tral energy distri bution (SED) 
of these IRAC counterparts, lYun et all (|2008l ) showed that 
SMGs as a population have characteristic red I RAC colors, 
simila r to dust obscured AGN as proposed by lLacv et al.l 
|2004h and IStern et alj (120051 ). These SMGs are systemat- 
ically offset from the color region associated with power- 
law AGN, however. Citing theoretical color tracks of dust 
obscured starbursts, Yun et al. advocated a dust-obscured 
young stellar population as the origin of the red IRAC color 
(see their Figs. 1 & 2). Objects with red IRAC color are rare 
(~ 1 arcmin -2 for [3.6] — [4.5] ^ 0.0) and distinct from the 
large number of foreground galaxies with characteristically 
blue IRAC colors. Both of these qualities can be success- 
fully exploited for distinguishing the SMG counterpart can- 
didates. Several color combinations are proposed by Yun et 
al., and we adopt here the simplest form, [3.6] — [4.5] ^ 0.0, 
since these two bands are the most sensitive and the most ro- 
bust among the four IRAC bands. In computing P-statistics 
for the IRAC counterparts, we thus use the number density 
of IRAC sources with [3.6] - [4.5] ^ 0.0, {S/N) 3 . 6 > 4, and 
S3. 6 > lAtJy- These are 36900 deg -2 and 31400 deg' 2 for the 
GOODS-S and SIMPLE IRAC catalogs, respectively. 

Examining the Spi tzer IRAC and MIPS properties of 
73 radio-selected SMGs. lHainline et al.l (|2009T ) reported that 
an IRAC color selection method similar to what we adopted 
is more successful in identifying correct counterparts than 
the IRAC c olor-m agnitude selection method described by 
I Pope et al.l (|2006l ). but they caution that the density of 
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sources meeting the lYun et alj (|2008T l color selection crite- 
ria is high enough to diminish the utility of this method. 
We adopt a more selective limit of [3.6] — [4.5] ^ 0.0, which 
is more effective in reducing the foreground confusion. In 
addition, we also employ a P— statistic analysis to give our 
method a more discriminative power. 



2.2 Counterpart Identification Results 

Finding charts for the 48 AzTEC/GS sources (GS2 is split 
into two) are show n in Fig. [1] in the order they appeared in 
IScott et ail l^Old ). Sources meeting the radio, MIPS 24 urn, 
and red IRAC color selection criteria are identified in each 
30" x 30" image centered on the AzTEC source position. A 
unique counterpart is easily identifiable in about 50% of the 
cases while two or more candidates are present in others, 
requiring a more systematic and quantitative analysis. 

Candidate radio and Spitzer counterparts and their 
computed P— statistics are given in Table [1] All "robust" 
counterparts satisfying P ^ 0.05 in any of the 3 bands are 
highlighted in bold-faced letters based on the analysis of the 
radio (Pi a), MIPS 24 urn (P 2 4 M ), or the IRAC color (P co ior) 
properti es. For 21 out of 48 cases, an 870 um LABOCA 
source (|Weiss et all [2009) is found within a 10" radius. 
Given the extremely low source density at 870 vim, the like- 
lihood of a chance coincidence is essentially zero. Therefore, 
we elevate the status of the 13 AzTEC/GS sources initially 
classified as only as a "tentative" identification based on the 
3 bands analysis to "robust" by folding in the astrometry of 
the 19" resolution 870 um LABOCA Survey data (see Ap- 
pendix A and Table [T]) - the remaining AzTEC-LABOCA 
sources are already classified as "robust" . We note that the 
low rate of coincidence between the LABOCA and AzTEC 
surveys (21/48 = 44%) can be largely accounted by the low 
S/N detections of sources in both surveys, although the pres- 
ence of high rcdshift (z > 3) so urces detected by AzTEC at 
1.1mm (e.g., lEales et al.ll2003t ) may play a role. Taking ad- 
vantage of the available rich multiwavelength database, we 
provide complete photometry for each source in Tabled 

A robust counterpart is identified for 27 (56%) out of 
48 AzTEC/GS sources using the P— statistic analysis com- 
bined with the LABOCA comparison. A total of 13, 8, and 
5 AzTEC sources have a robust counterpart based solely on 
the radio, MIPS 24 um, or IRAC color analysis, respectively. 
An additional 19, 14, and 18 have tentative identifications 
with 0.05 < P < 0.20, respectively. The robust radio and 
MIPS 24 vim identification rates of 13/48 (27%) and 8/48 
(17%) are consistent with other similar studies. For exam- 
ple, using similar depth Spitzer dat a and slightly dee per 
radio data in the GOODS-North field. IPope et al.l (|2006l ) re- 
ported robust identification rates of 21/35 (60%) and 6/35 
(17%) for the 1.4 GHz radio and MIPS 24 um data and ad- 
ditional 10 and 6 tentative identifications, respectively. For 
the SCUBA Half D egree Extragalactic Survey (SHADES), 
llvison et all (|2007l ) reported 56% and 54% robust identifi- 
cation using much shallower MIPS 24 urn and comparable 
depth radio data0The frequency of robust counterpart iden- 



tification rate using red IRAC color is similar to the MIPS 
24 urn identification rate, indicating that their respective 
candidate source density is comparable. 

Five AzTEC/GS sources (10%) have two robust coun- 
terparts each. This multiple robust candidate identification 
rate is similar to those found in the GOODS-North field 
JPope et al l 120061; IChapin et all 120091) and th e SHADES 
fields (|rvison et al.l [20071 : IClements et al.|[20oj ). This mul- 
tiple identification frequency is about 40 times higher than 
one would expect at random. A distinct possibility is that 
AzTEC c ounterpart sources ar e intrinsically strongly clus- 
tered (see I Williams et al.j|201ll . and references therein), and 
the P— statistic computes implicitly the likelihood that a 
particular candidate is either the AzTEC counterpart or a 
close companion. A strong clustering of SMGs is also ex- 
pected if they represent a rapid build-up of stellar mass for 



5 The MIPS 24 um robust detection rate by llvison et al.l j2007l l 
is significantly higher than ours or by IPope et alj d2006l') . despite 
their much shallower data, and this analysis may be in error. 



halo - see discussions bv lBlain et al.1 |2004). Similar SEDs 
and redshifts of the multiple candidate counterparts for 
AzTEC/GS19, AzTEC/GS31, and GN19 (|Pope et al.ll2006l ) 
offer further anecdotal evidence for the clustering explana- 
tion. IWang et al.l (|201ll ) has reported two examples where a 
single SMC is broken up into multiple discrete components 
when observed at high angular resolution with an interfer- 
ometer, further supporting the clustering sce nario. Based on 
the sim ulations of two large gas-rich galaxies. lHavward et al.l 
l|201ll ) have suggested that some fraction of SMGs may be 
such closely interacting pairs just prior to a merger, and such 
a scenario would certainly boost the pair fraction. Uncer- 
tainties in the parameters chosen for the simulations, such 
as the details of the progenitors and the microphysics of 
star formation and gas consumption, make the comparison 
with the observed statistics difficult. These new observations 
should serve as important observational constraints for fu- 
ture modeling studies. 

As shown in Table [2] only 22 out of 47 AzTEC/GS 
sources have an 870 um LABOCA counterpart in the pub- 
lished catalog bv lWeiss et al.l (2009). Utilizing the radio and 
MIPS 24 micron P-statistic and the IRAC 3.6 and 5.8 urn 
color- magnitude selection bv lPope et al.l (|200 6). Biggs ~et al.l 
(2011) identified 16 secure and 3 tentative counterparts 
among these 22 sources in common. In comparison, we iden- 
tify 16 robust counterparts based on the P-statistics alone, 
and all 19 individual candidates identified by Biggs et al. 
is either a robust (8) or tentative (11) counterpart in Ta- 
ble Q] The agreement between our results and theirs is very 
good mostly because both groups rely heavily on the radio 
continuum data for the counterpart identification. 



2.3 Counterpart Identification for SMGs in 
GOODS-North 

To improve the statistics of the subsequent analysis, we also 
apply the same counterpart identification methods to the 
AzTEC 1.1mm sources identified in the GOODS-North field 



Although the majority of radio-identified sources in Table [T] as 
well as by Pope et al. have a MIPS 24 urn counterpart, the MIPS 
24 urn source density is also much higher than the radio, leading 
to a greater chance-coincidence and thus a higher P— value in 
general. 
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Figure 1. Finding charts for the 48 AzTEC/GS sources. These false color images are 30" X 30" in size and produced using the IRAC 
3.6 urn (blue), 4.5 (i.m (green), and 8.0 \J.m (red) band images. Red circles mark the 1.4 GHz radio sources while yellow squares are MIPS 
24 u.m sources. For sources without a plausible radio or MIPS 24 u.m counterpart, IRAC sources with red IRAC color ([3.6] — [4.5] > 0) 
are identified using magenta squares. 
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Table 1. Radio and Spitzer identifications of AzTEC sources (procedure described in Section [JJ- The counterpart search radius Rs 
and the likelihood P values are described in detail in the text (§0, and robust counterparts are emphasized in boldface. Spectroscopic 
redshifts are given in the column labeled z spec (references for these measurements are given in Appendix lAt . 
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(a) Radio sources identified with IRAC positions priors; (b) MIPS 24 fim flux and P-statistic determined from the FIDEL catalog; (c) 
IRAC fluxes and P-statistic determined from SIMPLE catalog; (d) robust identification based on AzTEC+LABOCA analysis. 
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Table 1 - continued 



AzTEC 
ID 


Rs 
(") 


radio coordinate 
(J2000) 


Dist. 

(") 


Spitzer coordinate 
(J2000) 


Dist. 
(") 


PlA 
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[3.6] -[4.5] 
(mag) 
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-0.55 


0.735 


GS38b 
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8.0 


0.120 a 
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+0.23 
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J033302.71-275642.5 


8.6 
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0.132 
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GS41c 
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0.053 c 


+0.23 
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0.173 6 


0.208 c 


+0.15 




GS42a d 


6.9 






J033314.16-275612.0 


4.6 






0.146 c 


+0.04 




GS43a d 
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J033302.90-274432.9 


4.7 






0.157 c 


+0.23 




GS44a 


10.4 


J033240.84-273752.3 


9.3 


J033240.84-273752.6 


8.9 


0.151 


0.389 6 




-0.07 




GS45a 


12.2 


J033218.65-273743.3 


12.1 


J033218.58-273742.3 


12.0 


0.244 a 


0.620 6 




-0.32 




GS45b 








J033219. 09-273733. 5 


0.9 




0.004 6 


0.006 c 


+0.18 




GS45c 








J033219. 21-273731. 5 


1.9 




0.016 6 


0.026 c 


+0.21 




GS45d 








J033218.94-273730.0 


4.3 






0.129 c 


+0.20 




GS46a 


13.0 






J033157.27-275656.2 


6.2 




0.226 6 


0.255 c 


+0.17 




GS47a d 


12.2 


J033208.27-275814.0 


7.6 


J033208.23-275813.9 


7.8 


0.105 a 


0.280 6 


0.371 c 


+0.46 





l)Perera et al.ll2008l ). using the updat ed positions and pho - 
tometry for these so urces presented in iDownes et aH |201ll ). 
IChapin et all (2009) reported one or more robust counter- 
part to 21 out of 29 AzTEC sources and at least one ten- 
tative counterpart for the remainder. Our analysis, employ- 
ing slightly different criteria, identifies 16 robust counter- 
part sources and one or more tentative counterpart to all 
but two of the remaining sources. The agreement between 
Chapin et al. and our work is quite good, as 13 out of 16 ro- 
bust counterparts we identified were also identified as robust 
counterparts by Chapin et al. 



3 REDSHIFTS AND SPECTRAL ENERGY 
DISTRIBUTIONS 

Spectroscopic redshifts, z spec , are available for only a 
small subset (~ 30%) of candidate counterparts de- 
spite the extensive redshift surveys that have been 
conducted in the GOODS-South field over the years 



dLe Fevre et all 12004 


. Szokolv et all 2004: Mienoli et al.1 


120051; Vanzella et al 

\„ . , 1 — im 177- 


20051. 2006J; jNorris et al.l 


2006; 



Treister et al.ll2009l; IWuyts et al.ll20 09; Bale stra et aill2010l ; 
Silverman et al l |2010| ; Icasev et al.l [201 ll ). The" primary 



reason for this is that many of the candidate counterparts 
are extremely faint in the optical, with a median brightness 
of i ~ 25.4 among those listed in Table [2] (also see Fig. [6]). 
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Table 2. Photometry data listed in the same order as the id entifications in Table [T] All upper-limits are given at a significance of 
3-cr. De-boosted AzTEC 1.1mm flux densities are taken from iDownes et"al] l|201ll ). The LABOCA 870 |xm photometry comes from 
IWeiss et alj J2009T) . 



AzTEC 


1.4 GHz 


1.1 mm 


870 |xm 


24 [im 


8 |J.m 


5.8 H-m 


4.5 yim 


3.6 |xm 


i 


K 


ID 


(^Jy) 


(mjy) 


(mjy) 








t , T \ 

(^Jy) 


/ . T \ 


(mag) 


(mag) 


GSla 


32.0 ±6.3 


u 


7 +0.6 
' -0.7 


9.2 ± 1.2 


122.0 ±5.2 


28.2 ± 0.7 


20.0 ±0.6 


14.60 ±0.09 


10.39 ±0.06 


> 25.3 


> 22.9 


GS2.1a 


50.7 ±6.2 


O 


4 +0.7 
^-0.6 


9.1 ± 1.2 


148.0 ±4.1 


21.4 ±0.5 


14.0 ±0.4 


10.62 ±0.06 


7.50 ± 0.04 


26.1 


24.2 


GS2.1b 


44.1 ±6.2 






< 13.9 


< 1.3 


< 1.2 


< 0.19 


< 0.11 


> 27.5 


> 24.5 


GS2.1c 


< 18 








16.0 ±3.2 


6.9 ±0.5 


7.8 ±0.4 


7.73 ±0.06 


10.36 ±0.04 


22.1 


21.8 


GS2.2a 


< 18 


A 
1 


n +o.6 




62.9 ± 3.9 


16.2 ± 0.5 


15.6 ±0.4 


17.98 ±0.07 


22.28 ±0.04 


22.9 


21.1 


GS2.2b 


< 18 






83.2 ±4.2 


23.8 ± 0.5 


24.8 ±0.4 


20.46 ±0.07 


16.10 ±0.04 


26.4 


23.3 


GS2.2c 


< 18 








< 12.8 


14.3 ± 0.5 


19.9 ±0.4 


16.75 ±0.07 


15.30 ±0.04 


26.6 


21.3 


GS3a 


40.7 ±6.5 




O+0.6 
B -0.5 


8.8 ± 1.2 


49.2 ± 2.8 


13.0 ±0.5 


8.7±0.4 


6.17 ±0.06 


4.38 ±0.04 


> 25.3 


> 22.9 


GS3b 


< 19 






< 12.6 


3.1 ± 0.5 


4.2 ±0.4 


3.13 ±0.07 


2.76 ± 0.04 


> 25.3 


> 22.9 


GS4a 


25.4 ±6.5 


5 


1+0-6 
1 -0.6 


8.8 ± 1.2 


< 16.7 


18.6 ± 0.9 


13.6 ±0.7 


9.32 ±0.13 


7.21 ±0.06 


26.8 


> 24.5 


GS5a 


96.4 ±6.7 


4 


O+0.6 

°-0.7 


3.9 ± 1.4 


521.7 ± 10.9 


58.2 ± 0.9 


62.6 ±0.9 


73.70 ±0.18 


59.67 ±0.11 


23.1 


20.7 


GS5b 


111.7±6.7 






282.3 ± 7.7 


16.1 ± 0.9 


17.9 ±0.9 


14.12 ±0.18 


9.50 ±0.11 


25.7 


22.6 


GS6a 


31.0 ±6.3 


Q 
O 


fi +0.5 
b -0.6 


5.8 ± 1.4 


141.0 ±4.0 


27.4 ±0.5 


25.6 ±0.4 


17.70 ±0.06 


11.72 ±0.04 


28.3 


22.8 


GS6b 


< 18 






94.2 ± 3.7 


11.4 ±0.5 


12.1 ±0.4 


17.00 ±0.06 


21.09 ±0.04 


23.5 


21.1 


GS7a 


51.2 ±6.4 




O+0.6 
°-0.7 


9.1 ± 1.2 


103.0 ±9.3 


22.4 ±0.6 


17.7 ±0.6 


12.03 ±0.12 


7.89 ± 0.08 


> 25.3 


> 22.9 


GS7b 


< 19 






60.0 ±9.5 


5.0 ±0.9 


7.6 ±0.7 


9.33 ±0.23 


9.03 ± 0.09 


24.2 


> 22.9 


GS8a 


71.4 ±6.6 


■3 

o 


4+0-6 

^-0.6 


7.5 ± 1.2 


620.0 ±6.5 


42.9 ± 0.7 


62.6 ±0.6 


50.07 ±0.10 


36.85 ±0.06 


24.7 


21.6 


GS8b 


< 19 






164.0 ± 4.8 


27.5 ± 0.6 


20.6 ±0.5 


15.73 ±0.09 


13.44 ±0.05 


26.6 


> 22.9 


GS9a 


86.8 ±6.6 




D -0.6 




228.9 ± 10.3 


27.3 ± 0.9 


14.9 ±0.9 


12.37 ±0.17 


7.69 ±0.13 


25.3 


> 22.9 


GS9b 


< 19 






117.4 ±9.7 


3.8 ±0.9 


6.3 ±0.9 


7.27 ±0.17 


5.58 ±0.12 


> 27.5 


22.6 


GS9c 


< 19 








< 23.1 


8.7 ± 0.9 


8.0 ±0.9 


8.06 ±0.17 


5.99 ±0.13 


> 25.3 


> 22.9 


GSlOa 


89.3 ±6.4 


Q 

o 


O+0.7 
S -0.7 


7.6 ± 1.3 


214.0 ±8.4 


32.0 ±0.7 


40.5 ±0.9 


47.03 ±0.22 


41.21 ±0.15 


23.6 


21.3 


GSlOb 


< 19 








< 22.7 


4.8 ± 1.0 


6.9 ±1.2 


5.73 ±0.22 


5.52 ±0.15 


26.8 


> 22.9 


GSlla 


46.0 ±6.4 


Q 

o 


O+0.6 

°-0.6 




117.0 ±4.5 


32.5 ± 0.4 


23.8 ±0.3 


22.45 ±0.05 


22.89 ±0.04 


> 25.3 


> 22.9 


GS12a 


21.0 ± 6.5 


3 


i +0.6 
i -0.6 


5.1 ± 1.4 


31.6 ± 5.1 


12.2 ± 0.5 


7.2 ± 0.5 


3.89 ± 0.07 


3.54 ± 0.04 


25.2 


> 22.9 


GS13a 


22.8 ± 6.3 


3 


i +0.6 
i -0.6 




34.7 ± 3.6 


6.2 ± 0.5 


6.0 ± 0.4 


5.32 ± 0.07 


4.24 ± 0.04 


26.8 


22.5 


GS13b 


24.0 ± 6.3 






112.0 ± 3.8 


13.8 ± 0.5 


14.5 ± 0.4 


14.03 ± 0.07 


13.72 ± 0.04 


23.2 


22.0 


GS13c 


23.7 ± 6.3 








224.0 ± 3.7 


31.9 ± 0.5 


33.1 ± 0.4 


42.70 ± 0.07 


53.65 ± 0.05 


22.5 


20.2 


GS14a 


< 18 


2 


rv+0 6 
a -0.5 


3.3 ± 1.2 


< 12.8 


< 1.3 


< 1.2 


1.24 ± 0.06 


1.19 ± 0.04 


25.2 


23.9 


GS15a 


07 a _|_ c c 


3 


,+0 7 
q-ru. t 

a -0.8 


4.2 ± 1.4 


< 24.2 


16.0 ± 1.0 


13.2 ± 0.9 




fi 1 7 + n 1 2 


> 25.3 


> 22.9 


GS15b 


<^ 1 Q 






< 24.6 


5.8 ± 1.0 


5.7 ± 0.9 


^71 IS 


^ 47 + D 19 


> 25.3 


> 22.9 


GS15c 


< 19 








< 25.0 


< 2.9 


< 2.6 


1.12 ± 0.18 


1.12 ± 0.12 


> 25.3 


> 22.9 


GS16a 


30.7 ± 6.4 


2 


7+0.5 
'-0.6 




46.4 ± 3.3 


16.4 ± 0.5 


10.9 ± 0.4 


8.12 ± 0.07 


4.98 ± 0.04 


26.5 


25.8 


GS16b 


22.1 ± 6.4 






144.0 ± 3.5 


15.7 ± 0.5 


15.9 ± 0.4 


21.74 ± 0.07 


28.44 ± 0.04 


23.0 


20.6 


GS17a 


26.1 ± 6.3 


2 


n+0.6 
9 -0.6 




71.3 ± 9.6 


5.1 ± 0.4 


4.6 ± 0.3 


4.21 ± 0.05 


5.42 ± 0.03 


23.6 


22.6 


GS17b 


< 18 








61.9 ± 7.2 


21.2 ± 0.4 


26.3 ± 0.3 


20.19 ± 0.05 


16.94 ± 0.03 


26.0 


22.8 


GS17c 


< 18 








200.0 ± 5.3 


20.9 ± 0.4 


23.9 ± 0.3 


16.50 ± 0.05 


11.79 ± 0.03 


26.7 


22.7 


GS17d 


42.1 ± 6.2 








< 12.3 


11.3 ± 0.4 


12.0 ± 0.3 


9.07 ±0.05 


7.37 ±0.03 


25.2 


23.6 


GS17e 


37.9 ±6.2 








23.5 ±4.1 


9.7 ± 0.4 


7.1 ±0.3 


5.23 ±0.05 


3.88 ± 0.03 


27.6 


23.5 


GS18a 


25.1 ±6.3 


3 


1+0.6 
L -0.6 


6.4 ± 1.3 


84.6 ± 2.9 


13.1 ± 0.4 


11.7 ±0.3 


7.56 ±0.05 


4.91 ±0.03 


28.1 


> 24.5 


GS18b 


20.2 ±6.3 








126.0 ±4.0 


22.2 ± 0.4 


16.0 ±0.3 


10.85 ±0.05 


8.19 ±0.03 


24.5 


23.2 


GS18c 


< 18 








91.4 ±2.8 


11.2 ± 0.4 


12.6 ±0.3 


11.84 ±0.05 


16.97 ±0.03 


22.8 


20.8 


GS19a 


34.0 ±6.5 


2 


fi +0.5 
D -0.6 




317 ± 17 


18.9 ± 0.5 


25.8 ±0.4 


22.97 ±0.07 


18.74 ±0.04 


24.2 


21.7 


GS19b 


40.0 ±6.5 








149 ± 17 


15.9 ± 0.5 


22.7 ±0.4 


17.67 ±0.07 


12.72 ±0.04 


26.6 


21.7 


GS19c 


< 19 








105.1 ± 10.7 


9.1 ± 0.3 


12.3 ±0.3 


11.61 ±0.06 


9.97 ±0.04 


26.9 


21.6 


GS20a 


793 ± 99 


2 


7+0.6 
'-0.6 




4030.0 ±44.1 


3499.5 ± 0.6 


1113.8 ±0.6 


759.28 ±0.11 


1131.36 ±0.16 


16.1 


15.0 


GS21a 


43.6 ±6.3 


2 


7 +0.6 
'-0.7 




299.0 ±4.0 


16.6 ± 0.5 


24.0 ±0.4 


19.97 ±0.08 


15.98 ± 0.04 


25.3 


22.5 


GS21b 


< 18 








27.8 ± 2.6 


4.8 ±0.5 


6.9 ±0.4 


8.44 ±0.08 


7.51 ±0.04 


25.9 


> 24.5 


GS22a 


34.6 ±6.5 


2 


1+0.6 

1 -0.6 




290.0 ±4.4 


13.7 ±0.5 


15.3 ±0.4 


13.23 ±0.07 


10.02 ±0.04 


26.4 


21.5 


GS22b 


< 19 






< 11.4 


< 1.4 


2.7 ±0.4 


2.55 ±0.07 


2.26 ±0.04 


25.9 


22.8 


GS23a 


23.4 ±6.5 


2 


q+0.6 
°-0.6 




42.4 ±5.9 


19.7 ±0.5 


18.9 ±0.5 


12.19 ±0.07 


8.06 ±0.05 


24.5 


> 22.9 


GS23b 


35.2 ±6.5 






4.7 ± 1.4 


586.0 ±6.2 


36.4 ±0.5 


47.1 ±0.5 


46.86 ±0.07 


39.37 ±0.05 


24.0 


21.2 


GS24a 


18.5 ±6.1 


2 


,+0.6 
°-0.6 




90.8 ± 3.8 


16.4 ±0.4 


9.8 ±0.3 


6.61 ±0.05 


4.79 ± 0.03 


25.5 


24.1 


GS25a 


89.5 ±6.2 


1 


o+o : 
a -0.6 


5.9 ± 1.3 


140.0 ±3.6 


32.2 ± 0.5 


24.5 ±0.4 


18.81 ±0.06 


13.88 ±0.04 


23.9 


22.6 
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Table 2 — continued 



AzTEC 


1.4 GHz 


1.1 mm 


870 |im 


24 |im 


8 |im 


5.8 |im 


4.5 |im 


3.6 |im 


i 


K 


ID 


(MJy) 


(mjy) 


(mjy) 


(A«Jy) 




(A*Jy) 




(MJy) 


(mag) 


(mag) 


GS26a 


< 18 


9 9+0.5 
z ' z -0.6 




< 11.3 


2.4 ± 


5 


2.7 ±0.4 


3.64 ±0.07 


2.98 i 


04 


25.2 


24.6 


GS26b 


< 18 






< 12.1 


nnxn 
y.u ztz v 


K 
U 


7.1 ± 0.4 


6.71 ± 0.07 


J. OO HI u 


U 1 


24.9 


23.3 


GS26c 


< 18 






< 11.0 


11.5 ± 


5 


18.5 ± 0.4 


23.55 ± 0.07 


91 ceo + 


04 


25.9 


21.3 


GS27a 


23.6 ±6.5 


9 9+O.6 
z ' z -0.6 




171.0 ± 5.6 


25.3 ± 


5 


23.7 ±0.5 


16.57 ±0.07 


11.67 ± 


01 


24.8 


24.2 


GS28a 


< 18 


9 1+O.6 
A1 -0.S 




17.3 ±2.6 


10.8 ± 


5 


6.9 ±0.4 


4.33 ±0.06 


2.62 ± 


04 


26.5 


23.6 


GS28b 


< 18 




< 26.4 


s i 


n 
u 


< 1.0 


1.12 ± 0.06 


l.UJ ZLL U 




26.0 


24.6 


GS29a 


< 19 


9 q+0.6 

z -°-0.6 




54.4 ± 10.7 


1 e: q _1_ 1 
10. It 1 


A 

u 


27.4 ±0.9 


32.93 ±0.20 


4y .u / it u 


lo 


20.4 


19.8 


GS29b 


< 19 




47.4 ± 10.1 


O.Z zt 1 


U 


6.2 ± 1.0 


5.51 ± 0.19 


P. C\A _l_ n 
0.U4 it U 


lo 


23.8 


22.5 


GS30a 


37.2 ± 6.2 


i q+0.5 
i - 8 -0.6 




459.0 ±6.2 


9c; n _l_ n 
zO.y ± U 


tr 



35.4 ±0.4 


32.99 ±0.07 


97 K7 -L- n 
A 1 .01 it U 


n 1 


24.6 


21.0 


GS30b 


24.2 ± 6.2 






316.0 ±4.2 


97 1 " 
Z I .0 ± U 


r. 




33.7 ±0.4 


34.14 ±0.07 


zo.yo it U 


n 1 


26.3 


21.2 


GS30c 


< 18 






< 10.5 


9 A. -t- n 

Z.4 It u 


O 


2.7 ±0.4 


2.20 ±0.07 


1 RQ + 
i.uy it u 


1 
U 1 


25.7 


24.2 


GS31a 


25.1 ± 6.9 


9 9+0.7 
z - z -0.7 




427.0 ± 5.6 


31.5 ± 


7 


49.3 ±0.7 


51.71 ± 0.11 


fid qs + n 


07 


21.6 


19.8 


GS31b 


37.5 ± 6.9 






521.0 ±6.7 


f\n n -t- n 


e 
o 


75.4 ±0.7 


70.27 ±0.11 


Q7 AK _)_ n 

y / .40 it u 


07 


22.5 


19.6 


GS32a 


30.3 ±6.8 


9 q+0.8 

z -°-0.8 




371.1 ± 11.7 


32.4 ± 


9 


47.0 ±0.9 


45.47 ±0.17 


«7 09 + n 

U / . ZIZ u 


12 


21.6 


19.8 


GS32b 


34.8 ±6.9 




265.6 ± 11.9 


32.8 ± 


9 


35.8 ±0.9 


57.75 ± 0.17 


fid + 
iw.oy ziz u 


13 


22.4 


20.2 


GS32c 


46.3 ±6.9 






270.4 ± 12.2 


12.4 ± 


9 


12.2 ± 0.9 


11.06 ±0.17 


10.27 ± 


13 


25.4 


> 22.9 


GS33a 


28.6 ±6.2 


2 +0 - 5 




16.1 ±3.7 


5.1 ± 


5 


2.4 ±0.4 


1.70 ±0.07 


1.24 ± 


01 


26.9 


24.1 


GS34a 


33.0 ±6.3 


i 7 +0.5 
i -'-0.6 




173.0 ± 3.3 


32.9 ± 


5 


28.8 ±0.4 


35.68 ±0.07 


Q9 79 + n 
o£ . t zn u 


05 


23.4 


20.9 


GS34b 


< 19 






41.4 ±3.7 


14.3 ± 


5 


16.7 ±0.4 


20.34 ±0.07 


23.42 ± 


04 


23.5 


21.2 


GS34c 


< 18 






56.8 ±3.3 


5.8 ± 


5 


8.1 ± 0.4 


10.08 ±0.07 


8.83 ± 


04 


25.6 


22.3 


GS34d 


< 19 






< 10.5 


3.5 ± 


5 


4.7 ± 0.4 


4.54 ±0.07 


4.05 ± 


04 


24.3 


23.4 


GS34e 


< 19 






70.1 ± 3.2 


14.9 ± 


5 


17.2 ± 0.4 


19.86 ± 0.07 


17.31 ± 


04 


24.4 


22.1 


GS35a 


41.3 ± 6.7 


+0.6 

2 - 1 -0.6 




153.0 ± 3.8 


22.2 ± 


5 


19.8 ± 0.4 


14.72 ± 0.07 


10.51 ± 


04 


24.7 


23.1 


GS35b 


< 19 




< 12.9 


2.4 ± 


5 


3.1 ± 0.4 


2.76 ±0.07 


1.96 ± 


04 


25.3 


23.5 


GS36a 


< 19 


9 n + 0.6 




< 17.4 


7.6 ± 


7 


5.2 ± 0.6 


2.84 ±0.10 


1.52 ± 


06 


> 25.3 


> 22.9 


GS37a 


20.5 ± 6.4 


9 -i +0.8 
^•-"--0.8 




112.9 ± 7.0 


17.1 ± 1 





16.1 ± 1.0 


15.08 ±0.20 


13.95 ± 


13 


25.4 


> 22.9 


GS37b 


< 19 






< 21.5 


< 3 





3.6 ± 1.0 


2.89 ±0.19 


2.52 ± 


13 


25.9 


> 22.9 


GS38a 


220.0 ± 6.5 


t+0.6 
1 -'-0.6 




39.2 ± 2.6 


34.2 ± 


5 


58.1 ± 0.4 


67.61 ± 0.08 


112.41 ± 


06 


21.1 


18.8 


GS38b 


< 19 




184.0 ± 4.9 


15.8 ± 


5 


21.2 ± 0.4 


16.86 ± 0.07 


14.66 ± 


04 


> 25.3 


> 22.9 


GS39a 


37.2 ± 6.6 


1 -+0 7 
L -°-0.7 


3.8 ± 1.4 


< 26.9 


8.4 ± 1 





6.7 ± 0.9 


4.91 ± 0.19 


2.96 ± 


12 


> 25.3 


> 22.9 


GS39b 


< 19 






48.4 ±7.4 


5.9 ± 1 





8.5 ± 0.9 


9.17 ±0.18 


10.04 ±0 


12 


23.1 


22.3 


GS40a 


< 19 


1 q+0.6 
i ' 8 -0.7 




169.6 ±8.3 


7.9 ± 1 





13. 7± 1.1 


13.45 ± 0.20 


13.64 ±0 


13 


25.1 


22.0 


GS41a 


28.2 ± 7.0 


7 9+0.9 

'• z -i.o 




< 22.4 


4.9 ±0 


9 


6.5 ± 0.9 


4.35 ± 0.17 


3.52 ±0 


12 


> 25.3 


23.3 


GS41b 


223.0 ± 7.0 


12.0 ± 1.2 


62.0 ± 6.0 


4.7 ±2 


4 


6.6 ± 2.6 


5.35 ± 0.46 


4.06 ± 


34 


> 25.3 


> 22.9 


GS41c 


< 21 






< 23.9 


< 2 


7 


< 2.8 


1.71 ± 0.17 


1.39 ±0 


12 


> 25.3 


> 22.9 


GS41d 


< 20 






212.2 ±6.8 


17.1 ±0 


9 


19.9 ±0.9 


18.94 ±0.17 


16.55 ±0 


12 


24.3 


22.0 


GS42a 


< 21 


9 2+ 1 - 2 


14.5 ± 1.2 


< 26.6 


5.7 ±0 


9 


< 2.7 


2.50 ±0.18 


2.40 ±0 


13 


25.4 


> 22.9 


GS43a 


< 20 


6 1 +1 '° 


6.7 ± 1.3 


< 25.6 


4.0 ±0 


9 


< 2.7 


1.26 ±0.16 


1.02 ±0 


12 


> 25.3 


> 22.9 


GS44a 


43.5 ± 7.3 


o 9+O.8 
°- z -0.8 


t5.0±1.4 


143.1 ±9.3 


19.4 ± 1 





17.9 ± 1.0 


27.54 ±0.19 


29.39 ±0 


14 


23.7 


21.5 


GS45a 


33.7 ±6.9 


4-o±l:? 


8.1 ± 1.2 


446.3 ± 9.4 


91.9 ± 1 


1 


106.6 ± 1.1 


140.28 ± 0.21 


188.72 ±0 


15 


18.6 


17.8 


GS45b 


< 20 






64.4 ±9.3 


4.3 ± 1 


1 


5.9 ± 1.1 


5.02 ± 0.21 


4.24 ±0 


15 


> 25.3 


> 22.9 


GS45c 


< 20 






48.9 ±9.2 


7.3 ± 1 


1 


5.8 ± 1.1 


5.33 ±0.22 


4.41 ±0 


15 


26.4 


> 22.9 


GS45d 


< 20 






< 37.7 


3.8 ± 1 


1 


< 3.3 


2.11 ± 0.22 


1.75 ±0 


14 


> 25.3 


> 22.9 


GS46a 


< 21 


4.8tj'; 7 ' 


t4.8±1.4 


74.6 ±9.6 


6.7 ±0 


8 


9.0 ±0.9 


10.54 ±0.16 


9.01 ±0 


11 


25.1 


> 22.9 


GS47a 


43.2 ± 7.0 




7.3 ± 1.2 


72.3 ± 10.6 


9.7 ± 1 





9.6 ±0.9 


6.88 ±0.18 


4.51 ±0 


12 


> 25.3 


> 22.9 



t The AzTEC and LABOCA centroid positions are offset by a significant amount (> 10"). 



Robustly identified AzTEC counterparts are even fainter as 
discussed below. 

To learn more about the redshift distribution of these 
AzTEC sources and their nature, we rely on the extensive 
database of extremely deep, multi-wavelength broad-band 
photometry to analyze their SEDs using empirical and the- 
oretical models. We first examine the optical/near-IR pho- 
tometry data using standard methods for estimating photo- 
metric redshift ("photoz" hereafter), stellar mass (M„), and 
star formation rate [SFRjjv)- We also employ an indepen- 
dent analysis of the IR-to-radio SEDs to derive photometric 



redshift, IR luminosity, and dust-obscured star formation 
rate (SFR IR ). 

3.1 Optical/Near- IR SED Analysis 

The photometric redshift of each galaxy is computed 
by fitting the observed optical and near-IR spectral en- 
ergy distribution of the galaxies to stellar population 
synthesis models drawn from the PEGASE 2.Qj library 

B http://www2.iap.fr/pegase/ 
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Figure 2. Comparison of spectroscopic and photometric redshifts 
for the candidate AzTEC counterpart sources in the GOODS- 
South (Table [TJ field. Photo metric redshifts base d on the shal- 
lower MUSYC photometry llGawiser et alj |2006| ) are shown in 
empty symbols. The two dotted curves represent redshift uncer- 
tainties of Az/(l + z) = 0.036 (see the text). The distribution of 
Az/(1 + z) is shown as a function of z in the bottom panel. 

l|Fioc fc Rocca-Volmeran"gelll997l ) . The models are shifted in 
the redshift range of < z < 7 with a step size of Az = 0.01. 
For each galaxy, the weighted average photoz is derived as: 

fzP(z)dz 

Z pho to- jp {z)dz » W 

where P(z) is the probability distribution function of red- 
shift P(z) oc exp(— x 2 {z)). To evaluate the reliability of our 
photoz measurements, we compare the photozs with spec- 
troscopic redshifts (specz) of GOODS galaxies with reliable 
emission line redshifts in Figure [2] We find that the rela- 
tive error (defined as (z pho t — z spcc )/(l + z spcc )) has a zero 
mean (0.0005) and a very small deviation of 0.036 after 3<r 
clipping of the outliers. The fraction of outliers beyond 3<r 
is 9.9%. The means of the relative errors have no significant 
offset from zero at all redshift bins, especially for our inter- 
ested range of 1 < z < 5. The demonstrated accuracy of our 
photoz estimation is sufficient to justify a statistical study 
of the physical properties of our selected galaxies. Derived 
photometric redshifts of the AzTEC counterpart candidates 
are listed as z opt in Table [3] A blank entry notes that the 
optical counterpart is undetected or too faint. 

Physical properties (stellar mass, SFR, dust reddening) 
of the galaxies are measured by fitting the observed SEDs 
with the CB09 (Chariot & Bruzual, in prep.) theoretical 
stellar population synthesis models. The Salpeter IMF with 
a lower and an upper mass cutoff of O.IM© and lOOM© 
are adoptedQ The galaxy redshift is fixed to its z spec or 
z op t in this step. A regular grid of models spanning a wide 
range of E(B — V), star formation history (characterized 



7 Adopting the Chabrier or Krupa IMF will lower the der ived 
M* and SFR by a factor of 1.7 - see lSalimbeni et al.l i2009l') . 



by duration r and age ), and metallicity ar e examined. We 
apply the Calzetti law (|Calzetti et al.l [2000l for the internal 
du st extinction c orrection and follow the method described 
by iMadaul (|l995l ) to account for the IGM opacity. The \ 2 
value for each SED model fitting is computed as 

x — ^2 ' y^i 

where F oos>i , F mo d e i,i, and ct; are the observed flux, model 
flux, and observational uncertainty in the ith band. The nor- 
malization factor a is equal to stellar mass if F mo del,i is nor- 
malized to IMq in our pre-computed database. The model 
with the smallest \ 2 ls considered the best-fit model, and its 
parameters are used to compute the stellar mass and star 
formation rate SFRuv - see Table 

We estimate the systematic uncertainties in the derived 
quantities using simulations. We generate theoretical SED 
templates with different redshift, stellar mass, SFR, age, and 
dust extinction. In each band, we randomly draw a photo- 
metric uncertainty from the error distribution of all sources 
with the same magnitude from the GOODS parent photo- 
metric catalog and perturb the photometry of the template 
using a Gaussian random deviation with variance equal to 
the drawn photometric uncertainty. These steps are repeated 
100 times for each template in each band to generate mock 
SEDs. Our SED-fitting code is applied to these mock SEDs 
to derive the systematic uncertainties in the derived quan- 
tities. If redshift and IMF are known, a typical uncertainty 
in stellar mass from the SED-fitting is about 0.1-0.2 dex for 
all stellar masses. The uncertainty in SFR is about 0.1 dex 
for SFR > IOOMq yr _1 . If a photoz is used, a typical red- 
shift error of 8z/(l + z) = 0.05 translates to a 0.2 dex error 
in M* and SFR. Excluding the uncertainty in the IMF, the 
overall typical uncertainties in M* and SFR are 0.3 dex and 
0.5 dex, respectively. 

3.2 IR/mm/radio SED Analysis 

We derive an independent estimate of photometric red- 
shift, IR luminosity, and dust-obscured star formation rate 
(SFRm) by analyzing the observed IR/mm/radio part of 
the SED. First, a photometric redshift is derived using an 
up dated version of the photoz analysis method described 
bv lCarilli fc Yunl (|l999l ). Noting a remarkably tight correla- 
tion between radio and far-IR lumino sity for all star forming 
galaxies (see review bv lCondorJll992T ) and the rapid change 
in the observed flux density ratio between th e 850 urn band 
and the 20 cm radio continuum with redshift, ICarilli fc Yunl 
(1999) proposed this observed flux density ratio as a robust 
redshift indicator. The success of this method rests on the 
fact that the Rayleigh- Jeans (R-J) part of the dust spec- 
trum rises rapidly with frequency as S oc while the 
radio part of the spectrum falls as S oc u~ ' 75 , leading to 
more than two orders of magnitudes change in the observed 
flux density ratio between z = and z ^ 2. 

Subsequent anal ysis by iHughes et all (|2002f l and 
lAretxaea etail lj2003h have shown that incorporating ad- 
ditional photometric measurements in the far-IR to radio 
bands can improve the redshift estimate, but all of these 
methods are fundamentally limited by the intrinsic varia- 
tion in the SED, arising from variations in the nature of the 
energy source and geometry of dust distribution. To improve 
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Table 3. Redshifts and derived properties of the AzTEC GOODS-South sources. 
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AzTEC ID z« pcc z b opt logM, SFR UV z c MR SFR IR 

(M ) (M Q /yr) (Af /yr) 



GSla 




2.96 ± 0.45 


11.24 


94 


3.56+?-» 


439 


GS2.1a 




2.13 ±0.60 


9.81 


234 


3.20+^i! 


500 


GS3a 










3.091"? ii 




GS4a 




3.37 ± 0.25 


10.87 


75 


3 53+0.57 


416 


GS5a 


1.599 


1.66 ±0.60 


11.15 


1632 


2-03+1'! 

— U. ( o 


646 


GS6a 




2.47 ±0.65 


11.37 


37 


o 7O+0.60 
Z. O n no 
— 0.98 


220 


GS7a 


2.676 








2.56+O'g 


638 


GS8a 


2.252 


1.91 ±0.70 


11.04 


1553 


2.n+;;- 4 ,i 


466 


GS9a 




3.49 ±0.35 


11.01 


285 


i.98+;; : ^ 




GSlOa 


2.035 


1.58 ±0.55 


11.27 


1555 


2.03+° : « 


350 


GSlla 










2.5o+;; : ^ 




GS12a 


4.762 


4.55 ± 0.15 


11.28 


87 


o no + 0.70 

— 1.26 


803 


GS13a 




2.28 ±0.90 


10.20 


3172 


2.92+?'??, 
— 1.10 


250 


GS14a 


3.640 


3.50 ±0.30 


9.36 


459 


> 3.0 


600 


GS15a 




3.01 ±0.45 


10.00 


1024 


n qo+0.67 

— 1.13 


116 


GS16a 


1.719 


2.85 ± 0.60 


9.33 


114 


—0.95 


200 


GS17a 




1 m + n 1 n 

-L . V J -L _l_ \J . i-\J 


9.67 


3.4 


2 04+0.44 
z,a -1.08 




GS17b 




q 1 1 _i_ n 90 

O. ± 1 HI u.zu 


11 41 


11 


> 0.1 


oou 


GS18a 










J ' UU -1.14 




OS! Qa 




1 QQ _|_ n 


10.93 


99 


2 74 +0.52 
Z - ' 4 -1.04 


200 


GS20a 


0.037 


0.069 ± 0.038 


9.52 


339 


57+"' 1 ,! 
u u ' —0.41 


0.9 


GS21a 


1.910 


2.08 ± 0.65 


10.77 


295 


9 9 o+0.40 
° — 0.90 


322 


GS22a 


1.794 


2.42 ± 0.05 






9 qq+0.51 
z,c>9 -0.93 


204 












9 7 7+0.47 




GS23b 


2.277 


"i r 1 A i r\ or" 

1.64 ± 0.25 


11.35 


33 


9 oc+0.44 
z " 3o -0.90 


393 


GS24a 




1.94 ± 0.50 


9.37 


130 


„ ^ . 4-0 74 

3 - 04 Ii.io 


200 


GS25a 


2.292 


1.82 ±0.40 


10.11 


467 


. in on 

1 -°^-0.68 


300 


GS26a 




3.64 ±0.15 


9.61 


47 


> 2.6 


150 


GS27a 


2.577 


2.55 ±0.45 


10.02 


383 


2 47+O.57 
-0.99 


322 


GS28a 




3.29 ± 0.65 


9.95 


1140 


> 2.6 


177 


GS29a 


0.577 


0.48 ±0.15 


10.45 


4.0 


> 2.6 




GS29b 


2.340 


2.34 ±0.15 


10.36 


75 


> 2.6 


200 


GS30a 




1.51 ±0.25 


11.02 


6.0 


9 9Q+0-56 
z,zt, -0.88 


182 


GS31a 


1.843 


0.77 ±0.10 


11.96 


0.3 




222 


GS33a 




2.45 ±0.15 


10.08 


6.0 




200 


GS35a 




2.96 ±0.35 


10.74 


127 


Q7+0.53 
z ' z '-0.91 


508 


GS37c 










> 2.6 




GS39a 










2 01+ 51 
z - ui -0.85 




GS42a 




2.37 ±0.15 


9.56 


4.1 


> 4.5 




GS43a 










> 4.0 




GS45b 




2.07 ±0.25 


10.20 


12 


> 3.4 




GS45c 




2.91 ±0.35 


10.46 


156 


> 3.4 




GS46a 




1.67 ±0.10 


10.47 


11 


> 3.5 




GS47a 










2 5Q+0- 55 
z,Ja -0.95 





(a) Zspcc is a spectroscopic redshift. See §[A]for the individual references; (b) z" pt is a photometric redshift derived from the analysis of 
the optical/NIR SED - see §[XT] and (c) Zrm is a new photometric redshift derived from the radio and AzTEC 1100 nm photometry — 
see § 13.21 The listed uncertainty corresponds to a redshift range that includes 68% of acceptable fits. 



the accuracy of the derived redshift and its uncertainty, we 
adopted a three times larger set of SED templates, adding 
34 new sources with two or more photometry measurements 
in the R-J part (150^ < A < 1500/i) of the dust SED and 
at least one ra dio continuum measur ement, mostly from the 
new study bv lClements et ail |2010l ). We opted to use em- 
pirical templates of observed SEDs rather than a library 
of theoretical templates because there is growing evidence, 
such as the tightness of the radio-IR correlation, suggesting 
that nature favors a certain subset of SEDs. 



Another important addition is the use of Monte Carlo 
simulations to improve the handling of measurement errors, 
noise bias, and the template variations. A notable outcome 
is that the derived redshift uncertainties, listed in Table [S] 
are asymmetric about the mean millimeter-to-radio photo- 
metric redshift (zmr)- Citing the flat tening of the IR part 
of the SED with increasing redshift, ICarilli fc Yunl (|200dh 
have previously noted the asymmetry in the scatter of the 
"mean galaxy model", but in the opposite sense from the 
uncertainties in the derived zmr- This actually makes sense 
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Figure 3. Comparison of spectroscopic and millimeter-to-radio 
(MR) photometric redshifts for the AzTEC sources with a robust 
counterpart in the G OODS-South (Table [TJ and GOODS-North 
llChapin et alj|2009h fields. Three objects with clear evidence for 
AGN activity are identified as empty circles. The two dotted 
curves represent redshift uncertainties of Az/(1 + z) = 0.10. The 
distribution of Az/ '(1+z) is shown as a function of z in the bottom 
panel. 



since the asymmetry in the mean template and the uncer- 
tainty in the derived zmr should be in the opposite sense. 
This comparison thus shows that the common practice of 
quoting the redshift uncertainty using the (sub)mm-radio 
spectral index met hod based on the scatter in the Carilli fc 
, Aretxaga et al.ll2007l;lDannerbauer et al.l 



Yun template (e.g.. 



20081: IChapin et al 



20091 ) is in error. 



The derived zmr is in good agreement with z spec in 
most cases as shown in Figure[3]with Az/ (1 + z) ~ 0.1. This 
is not surprising given that the well-know radio-FIR correla- 
tion ap pears to hold among high-redshift IR-selected galax- 
ies (see llvison et al.ll2010l ; iLacki fc Thompson! l20ld ). This 
estimator may be more accurate for starburst-dominated 
SMGs since two outliers at z = 4.05 and z = 4.76 are known 
AGNs, and a similar "radio-excess" due to an AGN contribu- 
tion in the radio wavelength has been previously seen among 
other high redshift QSOs (|Yun et all l2002i ; lYun fc Carilu] 
2002). In some cases (e.g., GS2.2a, GS29a) the derived zmr 
is completely inconsistent with their z 3pec , primarily because 
of their radio non-detection. Given that the radio-FIR corre- 
lation holds for all other objects, a likely explanation is that 
their optical counterparts are mis-identified, as is expected 
to happen in a small fraction (?C 5%) of cases (see § 12. ip . 

Once the redshift of an AzTEC counterpart is deter- 
mined, its IR luminosity Lir can be estimated by adopt- 
ing an SED template most consistent with the observed 
FIR/mm/radio photometry data. For the ease of a di- 
rect comparison with optically derived SFRuv, we fix 
the redshift of each source to z spec or z opt . Then we use 
this Lir to de rive a dust-obsc ured star formation rate 
(SFRir - e.g., lKennicuttj|l99Sl ). We adopt a set of the- 
oretical SED templates for an ensemble of GMCs cen- 
trally illuminated by young star clusters bv lEfstathiou et al.l 



(2000) that are shown to provide a good fit to a wide 
range of IR-selected sources at di f ferent redshifts (e.g. , 
Efstathiou fc Rowan-Robinsonl 120031 ; IClements et~ai1 120081 : 
Rowan- Robinson! 120101 ) . Star formation rates derived from 
the IR SED fitting, SFRir, are computed directly from 
the best-fit model star formation history (also assuming the 
Salpeter IMF) and are summarized in the last column of 
Table [3] Note that z opt and zmr do not always agree well, 
particularly for the cases where the counterpart is not se- 
cure. Therefore SFRir is derived primarily for the securely 
identified AzTEC sources with a z spec or a well-determined 
z op t only. 



4 REDSHIFT DISTRIBUTION OF 
AZTEC-GOODS SOURCES 

4.1 Derived redshift distribution 

The deep multiwavelength data and the extensive spec- 
troscopic redshift surveys covering the two GOODS fields 
offer the best opportunity to identify millimetre- and 
submillimetre-bright galaxies and to construct the most 
complete redshift distribution yet. By utilizing the analysis 
of AzTEC sources in both GOODS fields discussed in § [3] 
we now have the opportunity to augment our understanding 
of the SMG redshift distribution with improved statistics. 

In Figures [4] and [5] we show the redshift distributions 
for the robust AzTEC counterparts in both GOODS fields. 
We use spectroscopic redshifts when available (22 robustly 
identified sources) and photometric redshifts otherwise, and 
we plot separately the distributions determined using z op t 
and zmr in both Figures. The redshift distributions based 



0.3 



GOODS (optical) 
GOODS (mm-radio)- 
Chapman + (2005) 




Redshift 

Figure 4. Redshift distributions of AzTEC sources in both 
GOODS fields based on the optical photoz [z pti 

N = 38; blank 

histogram) and the millimetric photoz (zmr, N = 74; hatched 
histogram) are compared with that of the 76 radio-identified 
SCUB A sources with spectroscopic redshifts by Chapman et al. 
j2005l . shaded histogram). The solid blue curve shown is a log- 
normal distribution as a function of (1 + 2) with a mean redshift 
of z u = 2.6 and <r = 0.2. 
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GOODS (optical) 
GOODS (mm-radio)- 

Aretxaga + (2007 




Redshift 

Figure 5. Redshift distribution of AzTEC sources in both 
GOODS fields based on the optical photoz (z op t; blank histogram) 
and the millimetric photoz {zmr\ hatched histogram) are com- 
pared with that of t he 64 robustly detecte d 850 |xm sources in 
SHADES survey by lAretxaga et al.l ||2007| . shaded histogram). 
Because of the shallow radio data in the Subaru/XMM Deep 
Field (SXDF), only the Lockman Hole sources are included for 
the Aretxaga et al. SHADES redshift distribution. The solid blue 
curve shown is a log-normal distribution as a function of (1 + z) 
with a mean redshift of z u = 2.6 and a = 0.2. 



on optical photoz (z op t) and millimetric photoz (zmr) are 
qualitatively in a good agreement with each other, while 
the two methods are each subject to potentially significant 
systematic uncertainties. These redshift distributions show 
that 80% of sources are at z > 2, with 60% just within the 
redshift range between 2.0 < z < 3.3. The relatively small 
number of robust counterparts with only a lower redshift 
limit (9 out of 74) assures that the median value of z m ed ~ 
2.6 is a robust estimate. In comparison, using a different 
redshift estimator and analyzing the pr operties of 29 AzTEC 
sources in the GOODS-North field onlv. lChapin et alj (120091 ') 
derived a median redshift of z = 2.7, in good agreement. 

The asymmetric redshift distribution of the AzTEC 
sources in the two GOODS fields shown in Figures [4] and 
[5] can be described reasonably well as a log- normal distribu- 
tion of the form 



/(*) 



[tn(l+£)-tnfl+£u)]^ 
5^2 



(1 + z)o\TtK 



The solid curve shown in both Figures corresponds to a log- 
normal distribution with jz m = 2.6 and a = 0.2 in ln(l + z). 
No attempt is made to derive the best fit values of Zp and 
a since some of the redshifts are only lower limits. Never- 
theless, these nominal parameters simultaneously describe 
the rapid drop-off on the low-z side and the long tail on the 
high-z side. 



4.2 Comparison with previous studies 

The redshift distribution of the SMG population (and thus 
their cosmic evolution) is still poorly understood. A com- 



parison of the redshift distribution derived from the AzTEC 
GOODS survey sources and those of of previous studies 
further illustrates this point. A comparison of the AzTEC 
sources in the two GOODS fields with that of the radio- 
selected SCUBA 850 urn sources by I Chapman et all (2005, 
Figure 2} gives an immediate impression that the two red- 
shift distributions are substantially different. In particular, 
the population of z ^ f.5 sources present in the Chap- 
man sample is missing in our sample while the AzTEC- 
GOODS sample shows a broader higher redshift tail. It is 
important to understand the underlying causes of this dif- 
ference since many studies have assumed that the redshift 
distribution derived by Chapman et al. is consistent with 



that of the SMG population as a whole (e.g. 


Coorav et al. 


2010 


;lDave et al. 2010; Narayanan et al.ll2010b 


;IVieira et al. 


2010 


: Almeida et al. 2011; Williams et al. 2011). As noted 



by Chapman et al. and further discussed below (§ 14. 3|) . these 
differences may be rooted in the use of radio selection for 
defining the Chapman sample. 

Another insightful comparison is made in Figure [5] 
by examining the redshift distribution of the GOODS 
AzTEC sources with the photometric redshifts of the 64 ro- 
bustly detected 850_um sources in the SHADES survey by 
lAretxaga et al.l l|2007h . The Aretxaga et al. redshift distri- 
bution is also largely missing the low-z population, and the 
agreement with our redshift distribution is better. If we take 
into account that 11 out 64 redshifts by Aretxaga et al. are 
only lower limits with zu m ?S 1-5, the redshift distribution 
for the AzTEC-GOODS and the SHADES Lockman Hole 
sources is in fact very good. 

A more quantitative comparison of the derived redshift 
distributio ns is made using the Ast ronomy Survival Analysis 
rASURV; lFeigelson fc Nelsol3ll985l ) package, which properly 
takes into account lower redshift limits|f| As summarized in 
Table[4] the only pairs of redshift distributions showing non- 
negligible probability of being drawn from the same parent 
sample are between zmr and z op t and between zmr and the 
photomet ric redshifts of the 850 urn sources in the SHADES 
survey by lAretxaga et all (|2007l ). On the other hand, both 
Gehan's Generalized Wilcoxon test and the Logrank test 
suggest that there is at most 3% probability that the Chap- 
man et al. redshift distribution is consistent with those of the 
GOODS AzTEC sources or the SHADES 850 urn sources 
in the Lockman Hole region as analyzed by Aretxaga et 
al. The Kaplan-Meier estimator gives the mean redshifts of 
< z M R >= 2.689 ± 0.112 and < z opt >= 2.516 ± 0.129. The 
inclusion of lower limits in the redshift for radio-undetected 



8 This analysis assumes that the censored data follow a similar 
distribution to that of the measured population. 



Table 5. The K-S test probability that two redshift distributions 
are drawn from the same distribution, derived using the cumula- 
tive distributions of the same samples as in Table [4] but excluding 
lower limits. 



Zopt 


zmr. Chapman05 


Aretxaga07 


Zopt 


0.845 0.250 


0.098 


ZMR 


0.014 


0.362 


Chapman05 




0.0014 
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Table 4. Comparison of mean redshifts and probabilities that two redshift distributions are drawn from the same parent distribution 
derived using the ASURV, as described in Section [4] The first row shows the number of redshift data points (Nd) and lower limits (Ni) 
used. The second row shows the Kaplan-Meier estimator for the mean and standard deviation of the redshift distribution. The last three 
rows show, first, the Gehan's Generalized Wilcoxon test probability for each pair, and second, the Logrank test probability. 





%opt 


ZMR 


Chapman05 


Arctxaga07 


(N d , JV,) 


(38, 0) 


(74, 9) 


(76, 0) 


(64, 23) 


< z > 


2.516 ±0.129 


2.689 ±0.112 


2.000 ±0.104 


2.695 ± 0.098 


%opt 




0.19 / 0.21 


0.030 / 0.026 


0.06 / 0.30 


Z MR 






0.0001 / 0.0001 


0.51 / 0.85 


Chapman05 








0.0000 / 0.0001 



sources appears to be the primary difference for these esti- 
mates (see below) . The mean redshift of the Aretxaga et al. 
sample is 2.695 ± 0.098, in a good agreement with that of 
< Zmr >• In contrast, the mean redshift of the Chapman et 
al. sample is 2.000 ± 0.104, significantly lower. 

The robustness of these analyses is further tested by 
conducting a K-S test for the same pairs of redshift distri- 
butions but excluding lower redshift limits. Again, as sum- 
marized in Table [5] the hypothesis that the Chapman et 
al. redshift distribution is identical to the zmr or Aretx- 
aga et al. redshift distribution can be rejected with better 
than 98% confidence while the same hypothesis for the other 
combinations of pairs cannot be rejected. Nevertheless, both 
the survival analysis and the K-S tests are giving us consis- 
tent results in that the Chapman redshift distribution is 
substantially different from the SMC redshift distribution 
derived by us using the GOODS AzTEC sources and that 
of the Lockman Hole 850 [im sources in SHADES survey by 
Aretxaga et al. 

4.3 Wavelength-dependent selection bias 

As discussed in some detail bv lChapman et all (|2005l ). their 
radio-selection for a spectroscopic redshift survey is intrin- 
sically biased toward low redshift galaxies and those with 
cold dust temperature. The observed radio flux density suf- 
fers from a strong positive k-correction, fading faster with 
increasing redshift than expected from the inverse square 
law. This means the majority of faint (Siaghz > 30-40 
^iJy) radio sources are at z ^ 1 and only extremely lu- 
minous (Pi. 4ghz ^ 10 24 W Hz" 1 ) radio sources are de- 
tecta ble at 2 > 2 (e.g. ISmolcic et al. 2008; IStrazzullo et al.l 
2010). Chapman et al. discussed this foreground confusion 
and removed about 10% of the sources from their sample, 
but the significant number of z ^ 1.5 sources remaining 
in their sample (but unseen in our sample and Aretxaga 
et al. sample) suggests that they under-estimated the fore- 
ground confusion. Chapman et al. also considered the ef- 
fects of rapidly declining sensitivity of the radio data with 
redshift, but they focused mainly on the dust temperature 
dependence on radio-IR correlation, assuming that the ma- 
jority of sources undetected in the radio bands are at the 
same redshifts as those detected. The significant high red- 
shift (z > 3) tail for the AzTEC and SCUBA sources seen 
in Figures [4] and [5] suggests that the highest redshift sources 
are missing in radio-selected samples like Chapman et al. 
because the depth of the available radio data is not suf- 
ficient to detect most SMGs at z > 3. Mapping the true 



redshift distributions of SMGs will require future complete 
spectroscopic redshift sur veys using instrumen ts such as the 
Redshift Search Receiver l|Erickson et al.ll2007l ) on the Large 
Millimeter Telescope or the Atacama Large Millimeter Ar- 
ray. 

These discussions of k-correction and dust temperature 
raise another important question as whether SMGs selected 
at 850 urn and 1100 urn (and as a natural extension at 250- 
500 |im by the Herschel SPIRE instrument) are systemat- 
ically different. Since the dust peak passes through these 
bands at different redshifts, some wavelength-dependent se- 
lection effects are expected. For example, identification of 
"drop-out" objects in these bands has been suggested as 
a m eans to identify the highest redshift SMG population 
(see iPope fc Chanl |2010| ). When we noted the low rate 
of coincidence between the LABOCA 870 [im sources and 
our AzTEC survey initially (see S 12.2) 1. this wavelength- 
dependent selection bias was one of the causes we explored, 
although we eventually concluded that the low completeness 
of the both surveys is the primary cause. The general agree- 
ment between the r edshift distribut i ons o f SCUBA 850 [im 
selected sources by lAretxaga et ail l|2007l ) and the AzTEC 
1100 um sources seen in Figure [5] suggests that the redshift 
distribution and SEDs of the SMG population is such that 
the sources identified at these two wavelengths are indeed 
similar. 



5 OPTICAL AND IR LUMINOSITY AND 
STAR FORMATION 

5.1 Not all SMGs are faint and red in the 
rest-frame UV and optical bands 

Although SMGs are a recently recognized class of galaxies, 
their relatively high density (0.1-0.5 arcmin -2 ) and high lu- 
minosity (Lir > 10 12_13 Lq) suggest that they represent a 
significant component of the general galaxy population and 
should play an important role in the over all galaxy evolutio n 
scenario. In the "down-sizing" scenario (|Cowie et al.lll99"6l ). 
more massive galaxies are thought to have been assembled 
earlier in cosmic history, presumably with a higher star for- 
mation rate (SFR). Luminous infrared galaxies (LIRGs) and 
ultraluminous infrared galaxies (ULIRGs) with SFR > 10- 
100 Mq yr _1 are the dominant co ntributor to the cos- 
mic star formation history at z ~ 1 (|Le Floc'h et al.ll2005l ; 
iMagnelli et al.l 2009, 201j]), and a significant contribution 
by SMGs with SFR of > 10 2_3 M© yr" 1 would represent a 
natural progression at z > 1. Massive galaxies with stellar 
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Figure 6. Measured iC-band versus i-band magnitudes of the 
robustly identified AzTEC counterpart sources (circles). Lines of 
constant colors with i — K = +0, +3, & +6 are shown. Circles 
with filled dots are AzTEC sources with Chandra X-ray detec- 
tion in the GOODS-North and GOODS-South fields while empty 
circles are the sources without X-ray detection. Crosses repre- 
sent the 17 " securely" identified counterparts to the SCUBA 8- 
mjy survey |rvison et all [2002) and the SCUBA Lens Survey 
llSmail et al .120021) sources. Asterisks are SMGs securely identified 
by interferom etric measurements at millimetre and submill i metre 
wavelengths jlono et ail l200rj ; lYouneer etail |2007| . l2008l 120091 : 
lHatsukade et al.ll20l6T l 



mass M* > 10 Mp are though t to be already in place by 
z ~ 2 (|van Dokkum et al.ll2008l ). Galaxies with even higher 
SFR might be found at higher redshifts. 

What would these SMGs look like in the rest-frame 
UV and optical bands? And how do they fit into the 
larger population of high redshift galaxies identified in 
those more traditional bands? Optical/UV size, morphol- 
ogy, and luminosity could provide an impor tant test for their 
origin as merger-driven sta rbursts (e.g., iNaravanan et a~ 



2010b) or large disk systems (Efstathiou & Rowa n- Robinson 



20031 : iKaviani et aU 12003 : lHavward et al.ll201ll) fueled by a 
high rate of cosm ological gas accretion (|Keres et a.1.1 120051 : 
iDave et alj|2010h . Early studies in the optical and near- 
IR suggested a diverse population of bright, modest red- 
shift (z < 1 ) and faint , high- r edshift (z ~ 2 ) galaxies, a s 
reported by (Lilly et all dl999h. iBarger ■ et all i 19991. l200d) 
llvison et all (|200fj| ). iFox et all (f2002D . and llvison et all 
(|2002l ). However, high resolution int erferometric imaging 
studies at millimeter wavelengths IIBertoldi et al.l 200C 



Fraver et all l200d : iLutz et al.l l200ll : iDannerbauer etal 



2002) have shown that the SMG counterparts are often 



undetected in the optical bands. A study of a large sam- 
ple of radio-ide n tified SMGs using deeper optical data by 
I Chapman et al.l (|200ll ) showed that their counterparts are 
indeed quite faint (/ > 25), and Chapman et al. concluded 
that dust obscuration makes these galaxies essentially invis- 
ible in the ultraviolet bands. This conclusion is n ot univer- 
sally accepted, however - see llvison et all l|2002l) . Interest- 



ingly, [ChaEman~etjjI] (|2005l ) targeted their own sample for 
spectroscopy using the Keck Telescopes and successfully ob- 
tained emission and absorption line redshifts for about 50% 
of their sample. 

A major motivation for this work is to clear up the 
confusion about the rest-frame optical and UV properties 
of SMGs by examining a robustly identified large sample 
with significantly improved statistics by taking advantage 
of the deep multiwavelength data available in the GOODS 
fields. In Fig. |S] we examine the rest-frame UV and optical 
properties of SMGs by plotting the measured i-band and 
K-band photometry of robustly identified AzTEC sources 
in the GOO DS-South field (th i s work) and the G OODS- 
No rth field dPerera et all 120081 : IChapin et aU l2009h - also 
see lPope et a.1.1 (l2006h. Sources i dentified by the S CUBA 8- 
mjy survey (llvison et al.l |2002|; ISmail et al.ll2004 ) and the 
SCUBA Lens Survey (|Smail et al.ll2002h are also shown for 
comparison. A remarkable result is that SMGs span a very 
broad range of brightness in both i- and K-banA (rest frame 
A = 240 nm & 630 nm at z = 2.5), spreading over 10 mag- 
nitudes, or a factor of 10 4 in flux density. The apparent 
brightness of the AzTEC GOODS sources by themselves 
span over 7 magnitudes with a median brightness of i = 25.3 
and K — 22.6 when the upper limits are taken into account. 
Although there is some overlap with the sources identified 
by the earlier SCUBA 8-mJy survey and the SCUBA Lens 
Survey, the AzTEC GOODS counterpart sources are sys- 
tematically fainter by ~ 3 magnitudes on average. Because 
earlier SMG identification studies relied on Tf-band data too 
shallow to detect the majority of the AzTEC counterparts in 
the GOODS-South field, this means earlier works may have 
missed or mis-identified the counterparts in many cases. The 
few but highly secure SMG counterparts identified recently 
using deeper optical and near-IR data and high resolution 
interferometric imaging in t he millimeter and submillimeter 
bands (shown as asterisks -llono et al.ll2006l:lYounger et all 
20071, l2008l . 120091 : lHatsukade et alj|20ld : IWang et alj|2007l . 



20111 ) have a brightness distribution more closely matching 



that of the AzTEC sources in both GOODS fields. 

We also deduce from the observed scatter in Fig. [S] 
that there is at least a factor of 10 variation in the in- 
trinsic rest frame optical luminosity among these SMGs. 
When viewed together with SCUBA-detected sources, these 
SMGs form a broad color track centered roughly around the 
i — K = +3 line (short-dashed line) , which is quite red com- 
pared with field galaxies. Some sources show a relatively flat 
color (i — K = +1) while there are others with extremely 
faint i-band upper limits and colors redder than i — K = +6. 
The scatter about the mean relation appears to increase at 
K > 22, but the source density is also higher at these fainter 
magnitudes. The full range of the scatter perpendicular to 
the mean relation is about 4 magnitudes in color. This large 
spread in color substantiates the earlier sugges tion that op- 
tical proper ties of SMGs are quite diverse (e.g.. llvison et all 
2000, 2002). However it cannot fully account for the up to 
~ 10 magnitude spread in their apparent brightness as a 
population. In other words, there is an additional factor of 
> 100 variation (or > 5 mag) in the apparent brightness 
of SMGs on top of the apparent differences in color, which 
may be due to variations in SED and extinction. Given their 
extreme luminosity, the SMG phase likely represents a brief, 
special moment during the rapid mass build-up phase (e.g., 
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iNaravanan et al.ll2010bl ; lHavward et al.ll201lh . If the major- 
ity of SMGs are seen in the redshift range between 1.5 and 
4.0 (see § H}, then the spread in the luminosity distance can 
account for about a factor of 10 in the apparent brightness 
variation. Therefore, the remaining factor of > 10 scatter in 
apparent brightness has to be accounted for by an intrinsic 
scatter in the rest-frame optical band luminosity. 

In the broader context of understanding star forming 
galaxies in the early universe, some of the AzTEC sources 
are bright in the rest-frame UV and optical bands and are 
already identified as star forming galaxies by past surveys. 
For example, about 30% are bright enough in th e optical and 
NIR bands to be classified as "BzK" galaxies (Daddi et al.l 
l2004h using the existing photometry - see discussions in Ap- 
pendix. Some of the z ~ 4 SM Gs have also been identified 
as "Lyman break galaxies" (see lCapak et al]|2008h . On the 
other hand, the majority of the AzTEC GOODS sources are 
too faint and red to have been identified in previous surveys 
of star forming galaxies and are likely to be entirely missed 
in the current accounting of the cosmic star formation his- 
tory. Future millimetre wavelength surveys with higher an- 
gular resolution are needed to probe deeper into the lower 
flux density (and lower luminosity) regime in order to bridge 
these populations and obtain a complete census of star form- 
ing galaxies. 



5.2 SMGs as massive galaxies in a phase of rapid 
stellar mass build-up 

5.2.1 Stellar luminosity of SMGs 

One constant in the high resolution interferometric millime- 
tre and submillimetre observations of SMGs is the presence 
of a Spitzer IRAC counterpart in the 3.6-8.0 urn bands, and 




Figure 7. Spitzer IRAC 3.6 um band "Hubble diagram". IRAC 
3.6 um band brightness (m3.6 M ) of the robust AzTEC counter- 
parts is plotted as a function of redshift. Symbols are the same 
as in Fig. [6] Small dots on the left panel are K— band selected 
star forming galaxies with known spectroscopic redshifts in both 
GOODS fields. The apparent brightness of a 40 Myr starburst 
population with a stellar mass of 10 9 and 1O 1O M0 is shown in 
dotted lines, while the same population after passively evolving 
for 500 Myr is shown in solid lines. An exponentially decaying 
starburst history with a 20 Myr e-folding time and solar metal- 
licity is assumed for the models. 



this is one of the key features we employ to identify AzTEC 
counterparts (see §[2|. A comparison of the apparent bright- 
ness in the 3.6 um band (rest-frame optical or near-IR) for 
the AzTEC counterparts and other K-band selected z ~ 2 
starforming galaxies with known spectroscopic redshifts in 
both GOODS fields is shown in Fig. [7^. While there is some 
overlap between these two populations, the AzTEC coun- 
terparts are systematically brighter by > 1 magnitude on 
average and represent the most luminous galaxies at these 
redshifts. This intrinsically high luminosity in the rest-frame 
optical and near-IR bands is clearly an important reason 
why these AzTEC sources are so readily detected by Spitzer. 

Since AzTEC detection requires highly efficient dust- 
processing of the UV radiation from young stars, the high in- 
ferred luminosity in the rest-frame optical/NIR bands seems 
surprising. After all, we just established in the previous sec- 
tion (§ 15. ip clear evidence for severe attenuation of UV ra- 
diation among many of these objects. A natural explanation 
for this apparent puzzle is found in the studies of the local 
ULIRG population. An imaging study of local ULIRG s in 
the near- and far-UV bands bv lGoldader et al] (|2002h has 
shown that activity traced in the UV bands is distributed 
over kiloparsec scales and is heavily obscured, particularly 
in the regions of the most intense starburst activities. When 
observed at redshifts of z ~ 2 to z ~ 4, these ULIRGs are ex- 
pected to be extremely red and faint in the observed optical 
and NIR bands (R-K = 4-6, K « 21-24), similar to the ob- 
served valu es for the AzTEC s ources shown in Fig. [6] At the 
same time, [Ch en et al.l ([2010) have shown that the stellar 
hosts of local ULIRGs are also extremely blue in rest-frame 
optical bands and are on average ~ 1 magnitude brighter 
than the field star-forming population, owing to distributed 
star formation activity and the high intrinsic luminosity of 
young stellar clusters. A H ubble Space Tele s cope NICMOS 
and ACS imaging study bv lSwinbank et all (|2010t ) has also 
found evidence for ongoing mergers and structured dust ob- 
scuration among z ~ 2 SMGs, further supporting the paral- 
lel in the observed source luminosity and structure between 
the local ULIRGs and high-z SMC population. 

One cannot automatically conclude from their large ob- 
served luminosity that these SMGs are also the most mas- 
sive galaxies at their observed epochs if their luminosity is 
powered by a large population of widely distributed young 
stellar clusters. The apparent 3.6 urn brightness of a sin- 
gle stellar population starburst model with a total accumu- 
lated stellar mass (M») of 1O 9 M and 10 lo Af Q is shown in 
Fig-El 3 for t wo different scenarios: 40 Myr (dotted lines) and 
500 Myr (solid lines) after the initial burst. The observed 
brightnesses of AzTEC sources are well-bounded by the 40 
Myr old starburst models with stellar masses of 10 9 M Q and 
10 10 M o , which are about 10 times smaller than the stellar 
masses derived for the K— band selected star fo rming galax- 
ies sh own in comparison (M» = 10 10 - U M Q , [Daddi et al.l 
2007). However, after just 500 Myr of passive evolution, the 
same starburst systems fade by ~ 2 magnitudes at z ~ 2-4, 
bringing them back to a better agreement with the mass esti- 
mates for the if— band selected star-forming galaxies. Alter- 
natively, accounting for the observed 3.6 um band brightness 
of the AzTEC sources assuming a maturing stellar popula- 
tion would require stellar masses well in excess of lO n M . 
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Figure 8. The stellar masses of AzTEC counterparts derived 
from optical and NIR photometry as a function of redshift. A 
histogram of stellar masses is shown on the right panel. Sources 
with a spectroscopic redshift are shown as solid symbols while the 
ones with a photometric redshift arc shown as empty circles, and 
they represent similar ranges of stellar mass. Typical uncertainties 
for the M„ and photo-z estimates are shown on the bottom left 
corner. Those detected by the Chandra in the X-ray are identified 
with a larger circle. Only GOODS TFIT catalog sources with a 
proper stellar mass estimate are included. 



5.2.2 Stellar mass and star formation rate of SMGs 

We can get a better handle on the stellar mass by modeling 
the observed rest-frame UV and optical SED as discussed 
in § 13.11 The derived stellar masses from the modeling of 
the UV-optical SED, shown in Figure [8] range between 10 9 
to 1O 12 M0. The majority of the derived stellar masses are 
between 1 and 30 times 10 10 M©, similar to those of the 
K — band selected massive star-for ming galaxies at th e same 
redshift, such as those discussed in lDaddi et alJ (|2007h . Sim- 
ilarly large stellar masses were found previously for SMGs 



(e.g-.lDve et alJl200StlTacconi et alJl2008l;lDaddi et al.ll2009l ; 
lHainline et al.ll2009l ; iMichalowski et al.ll2010T ). and they are 
consistent with the idea that these submillimetre-bright 
galaxies are associated with the peak of the stellar mass 
function at z = 2-3. 

Six out of 18 AzTEC counterparts with M* ^ 10 11 A/ 
are also X-ray sources detected in deep Chandra surveys. 
The frequency of the Chandra detected sources is about the 
same for the lower stellar mass galaxies (seven out of 25), 
and there is little evidence for any dependence on stellar 
mass. The low number of z ^ 3 sources detected by the 2-4 
Ms Chandra surveys may reflect the limiting sensitivity of 
the X-ray data, and the observed X-ray fraction is a lower 
limit. Given the poor statistics, it is difficult to conclude 
whether the presence of a luminous AGN is influencing the 
modeling of the rest-frame optical SEDs. 

The rest-frame optical SED modeling also yields a UV- 
derived star formation rate (SFRuv). The derived SFRuv 
for the AzTEC-GOODS sources cover a broad range: 1 — 
2000M© yr" 1 - see the left panel of Figure |U A sur- 



prising result is that the derived SFRuv is quite high, 
^ 100 - IOOOA/q yr" 1 for about 50% of the cases. The ob- 
served SFRuv distribution is also nearly independent of 
stellar mass. The SFRuv distribution broadly overlaps the 
observed SF R — M, relation for t he if-band selected galax- 
ies studied bv lDaddi et alJ |2007l . dashed line), but there is 
little evidence that these AzTEC sources follow the same 
SFR - M» relation. The AzTEC-GOODS sources also do 
not follow the SFR — M* relation predict ed by the z = 2.5 
model SMGs fueled by cold flow accretion (|Dave et al.ll2010l 
solid line). Some of the galaxies with the highest SFRuv 
are detected in the X-ray by Chandra, raising the possibility 
that the UV light from the central AGN might contribute 
to these high values. However, not all Chandra detected 
sources are associated with a high SFRuv , and neither the 
SFRuv / SFRir ratio nor the specific star formation rate 
discussed below offer any clear evidence to support this idea. 

For comparison, IR-derived star formation rates 
(SFRir) determined from modeling the IR SED are shown 
in the middle panel of Figure The SFRir is uniformly 
high, ^ 100 — IOOOM0 yr -1 , with a much smaller scat- 
ter and completely independent of stellar mass. This is ex- 
pected since these confusion-limited AzTEC surveys prefer- 
entially select sources with intrinsically large Lir. We note 
that the 1.1mm selection does not guarantee a high Lir or 
SFRir if cold dust (T<j = 10 — 20 K) emission dominates 
the millimetre spectrum. On the other hand, our SED mod- 
eling does not find any cold dust dominated sources with 
Lir < lO 12 I/0. Since the IR luminosity accounts for the to- 
tal amount of dust-processed UV radiation, a comparison of 
SFRir with SFRuv should offer a crude measure of the ge- 
ometry between the young stars and the obscuring dust. The 
large SFRuv derived for a large fraction of AzTEC-GOODS 
sources is particularly interesting in this regard, and this 
result may indicate that star formation activities and dust 
distribution in these SMGs are not as co ncentrated as in the 
local ULIRGs, where Lir » L bo i (see ISanders fc Mirabel 
Il996t) . 

To explore the relationship between SFRuv and 



1000 
100 
s 10 



K 
fc 

in 



fc 

m 



- a 



H 



1 

0.1 
0.01 
0.001 



(b 



□ 

I B D n 





a 3 
z 



.- 



1 1 1 
10 



1 n 



J L 

1 1 



I I 

12 



M. (MJ 



Figure 10. The ratio of SFR, derived from the rest-frame UV 
(SFRuv) an d the rest-frame IR (SFRir) as a function of (a) 
redshift and (b) stellar mass (M„). Galaxies with spec-z (photo- 
z) are shown in filled (empty) symbols, and those detected by 
the Chandra in the X-ray are identified with a larger square. The 
dashed line marks the SFRuv = SFRir relation. 
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Figure 9. Star formation rates for AzTEC GOODS sources estimated from the rest-frame UV (SFRuy) and IR (SFRir) as a function 
of stellar mass. Open and filled symbols represent photometric and spectroscopic redshifts, and those detected by the Chandra in the 
X-ray are identified with a larger symbol. A typical overall uncertainty for an object with photoz is shown on the bottom right corner 
of the left panel. The mean SFR - AT, relation for the z = 2.5 model SMGs fueled by cold flow jPave et al.l|20ld) is shown by a solid 
line. Long-dashed line is the observed mean SFR — M t relation for the i\-band selected galaxies iDaddi et al.ll2007l) . The SMGs do not 
seem to follow either trends. 



SFRir further, their ratios are plotted as a function of red- 
shift and M» in Figure [TO] This ratio varies widely from one 
source to another, spanning over 5 decades in total range, 
and it is independent of z and M*. Finding a large number 
of sources with SFRuv / SFRir > 1 is particularly puz- 
zling for these galaxies with a large stellar mass. A mis- 
identification of the counterparts is also a plausible expla- 
nation, but the observed distribution would require the fail- 
ure of counterpart identification in a large fraction of cases. 
Either an under-estimate of SFRir or an over-estimate of 
SFRuv (and possibly both) can provide an explanation, 
as the estimates of both SFR and M* ar e subject to sig- 
nifica nt systematic uncertainties (e.g., see iMaraston et al.l 
l2010l ). If these galaxies represent young galaxies seen during 
their rapid mass build-up phase (see below), then the well- 
known mass-metallicity relation (|Tremonti et al.ll2004l ) and 
the attenuation of UV light in the local universe may not be 
directly applicable. The presence of an X-ray source detected 
by Chandra is not correlated with the SFRuv / SFRir ra- 
tio, and the presence of an X-ray emitting AGN does not 
seem to contribute directly to the derived SFRuv in most 
cases. 



5.2.3 Specific star formation rate and mass build-up 
history 

Charting the star formation and stellar mass build-up his- 
tory is one of the most powerful tests for galaxy evolu- 
tion theories. For example, the emergence of red sequence 
galaxies that are massive and passively evolving around 
z ~ 1 and their increase in number with time are widely 
cited as important observational constraints that require 
additiona l complexities such as AGN feedback an d "dry" 
mergers (|Bell et alJliooi . 120071 ; iFaber et alj|2007l ). Statis- 
tical studies such as the Sloan Digital Sky Survey have 
shown that the bulk of stars now in massive galaxies formed 
at earlier epochs than stars in lower mass galaxies (e.g., 



iKauffmann et al.l l2003r i. suggesting a strong link between 
galaxy mass and star formation history. A particularly useful 
quantity to examine in this regard is the specific star forma- 
tion rate (SFR per unit stellar mass; SSFR = SFR/AQ. A 
systematic dependence of SSFR on galaxy mass and a rapid 
increase of SSFR w ith redshift have been established by sev - 
eral recent studies (|Zheng et al.ll2007l ; TPamen et ai1l2009h . 
Given the large stellar masses (M* = [1-30] x 1O 1O M ) and 
SFRs (> 10 2_3 M Q yr _1 ) for these z = 2-4 SMGs, examin- 
ing their SSFR in the context of the observed trends with 
stellar mass and redshift may provide a valuable new insight 
into the physical mechanisms driving the SMG phenomenon 
and massive galaxy formation. 

The computed SSFRs for the AzTEC-GOODS sources 
are shown in Figure [TT] An immediately noticeable trend is 
that the derived SSFRs are uniformly quite high, SSFR ~ 
1-100 Gyr -1 . Among the optically selected samples, galax- 
ies with SSFR > 10 — 100 Gyr -1 are generally associated 
lower stellar mass (Af* ^ 1O 1O M0) galaxies undergoing a 
starburst episode. More massive galaxies in the local uni- 
verse ar e associated with 1 -2 orders of magnitudes lower 
SSFR dBauer et all 120051 : iFeulner et ail 120051 ; lErb et all 
l200rj:lElbaz et al.ll2007f ) . The same trend also holds at higher 
redshifts as the if— band selected star formin g galaxies at 
z ~ 2 have on average SSFR « 1 Gyr -1 (Da ddi et al.l 
2007), overlapping only at the botto m range of the SSFR 
associated with the AzTEC sources. iPannella et~al] (|2009h 
have also reported an average SSFR w 5 Gyr" 1 and 
SFR » 100M Q yr" 1 for their z ~ 2 radio-identified star 
forming galaxies with an av erage M, = 3x 1O 1O M at 2 ~ 2 



rorming galaxies with an avei 
(also see lDunne et al.ll2009l ). 



The SSFRs for the AzTEC-GOODS sources are sig- 
nificantly higher than those of similar stellar mass galaxies 
in the local Universe, and they appear to follow the same 
broad trend of rapidly increasing SSFR with redshift. The 
dotted line shown in Figure [T T] is the stellar mass-d ependent 
SSFR evolution mapped by Damen et al.l (|2009h for mas- 
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Figure 11. Specific star formation rate of AzTEC counterpart 
candidates derived from the (a) rest-fra me UV and optic a l SEP 
fitting and (b) far-IR SED fitting using lEfstathiou et al.l feOOCh 
dusty starburst SED templates. Solid and empty symbols rep- 
resent the sources with spectroscopic and photometric redshifts, 
respectively, and those detected by Chandra in the X-ray are iden- 
tified with a larger symbol. A typical uncertainty is shown for one 
of the photo-z sources at the bottom of the left panel. The aster- 
isks represent the radio-derived SSF Rs for star forming ga laxies 
with M* ~ 3 X 10 10 M Q derived bv IPannella et alj (|2009h . The 
solid line represents the inverse of the Hubble time, and sources 
above this line are in a starburst mode. The dotted line connect- 
ing the crosses represents the average SSFR for m assive galaxies 
with M, 3= W 11 M Q found bv lDamen et alj J2009t) . 



sive galaxies with M* 10 1 Mq, and it shows a rapid 
rise as SSFR oc (1 + z) 5 between z = and z = 2. 
The radio-derived SSFRs fo r star forming galax ies with 
M,~3x 10 10 Af© derived bv lPannella et al'l (|2009l ). shown 
in asterisks, extend this rapidly rising trend to z ~ 2.5. The 
AzTEC-GOODS sources extend this rise in SSFR further 
to z ~ 4, although there is significant scatter. As noted 
by Damen et al. and others, there is a mass-dependence on 
the SSFR evolution, and the spread in M* for the AzTEC 
sources likely contributes to some of the observed scatter. 

The SSFRs derived for the AzTEC-GOODS sources 
provide the strongest evidence yet that SMGs are seen dur- 
ing the brief phase of rapid stellar mass build-up. The solid 
lines in Figure [TT] represent the inverse of the Hubble time, 
1/tH' only galaxies with SSFRs above this line have suffi- 
ciently high SFRs to build up their current stellar masses 
within the Hubble time at their respective redshifts. Consid- 
ering the SSFRs derived from fitting the FIR SEDs (right 
panel of Figure [TTj) , all of the AzTEC GOODS SMGs are 
located at or above this critical line. We cannot exclude the 
possibility that these SMGs are rejuvenated galaxies, under- 
going another episode of extreme luminosity, but it would 
require an even earlier episode of rapid stellar mass build- 
up. Citing extremely high SFR and similar density, previ- 
ous studies have made plausible arguments for identifying 
SMGs as progenitors of pr esent day massive elliptical galax- 
ies (e.g.. lBlain et al.ll2004h . Our new analysis of the SSFR 
allows us to demonstrate quantitatively that these SMGs 
are seen during a phase of rapid stellar mass build-up. 

The absence of AzTEC sources with SSFR below the 
1 /tu line in the right panel of Figure [TT] is primarily the 
result of AzTEC survey depth - in fact, all existing con- 
fusion limited surveys carried out with, e.g., AzTEC and 



Herschel, probe only the brightest end of the luminosity 
function. Much of the cosmic IR background (CIRB) is ex- 
pected to arise from fainter sources below the confusion 
limit, and their number counts can offer an important con- 
straint to the evolution model for SMGs jGranato et all 
12004 iBaugh et al1l2005l ; lRowan-Robinsonll2009T l. The loca- 
tion of the 1/tH line in Figure [TT1 leaves a fairly limited pa- 
rameter space for a lower luminosity dust-obscured starburst 
population that can contribute significantly to the CIRB - 
e.g., SFR ^ 20-100 M Q yr" 1 for a (1-5) xlO lo M galaxy 
at 2 = 2. A a > 2 galaxy with M* ^ 10 10_11 M© can still ap- 
pear with a SSFR below this 1/tn line, but the presence of a 
large population of such galaxies would have an important 
consequence in that the formation epoch of those massive 
galaxies has to be pushed to a much earlier time. The de- 
creasing tn with redshift also requires an even larger SFR 
with increasing z, and in turn the submm/mm-detected frac- 
tion of galaxies with a stellar mass M, ^ 10 10_11 M Q has to 
rise with increasing redshift. The high detection rate of op- 
tically s elected z ^ 4 QSOs in the submm/mm co ntinuum 
(~ 30%. ICarilli et al.ll200ll ; IWang et alj|2007ll2008h . despite 
the selection bias against obscured systems, appears to be 
in line with this expectation. 

The SSFR analysis of the AzTEC-GOODS sources also 
suggests an intriguing idea that there may be two classes of 
SMGs, possibly driven by two different modes of star for- 
mation or observed at two different phases. While the M* 
for these SMGs span over a factor of 30 (see Fig. [8}, the 
SFRir and SFRuv + SFRir show no depend ence on M* 
in Fig ured This is in contrast to the finding bv lDave et al.l 
l|201lf ). where a tight M*-SFR relation is a generic out- 
come of all of their cosmological hydrodynamic simulations 
incorporating galactic outflows. One way to interpret our 
observational results summarized in Figure [9] is that only 
SMGs with M, > 10 10 ' 5 have properties similar to the ob- 
jects modeled by Dave et al., following the M* -SFR rela- 
tion predicted for SMGs (solid line in Figure [9] iDave et al.l 
l2010h . AzTEC-GOODS sources with M» < 10 10 ' 5 may fol- 
low a different, currently unknown process that leads to 
10 times larger SSFR. This is somewhat analogous to the 
mass-dependent SSFR seen among galaxies in the local uni- 
verse ("red" and "blue sequence") with a similar range of 
M t , although with 1-2 orders of magnitudes lower SSFRs. 
A division by mass certainly seems somewhat arbitrary. On 
the other hand, these SMGs appear to show a sign of group- 
ing by mass in the middle and right panel of Figure [TT] A 
hint of bimodality in the histogram of M* is also seen in 
Figure [S] Future LMT and ALMA surveys of much larger 
samples with higher angular resolution and spectroscopic 
redshifts should provide a definitive test on this intriguing 
possibility. 



6 AGN AND STARBURST ACTIVITIES 

Determining the source of enormous luminosity associated 
with SMGs is an outstanding problem that has important 
implications on understanding the mass assembly history 
of galaxies. To probe the nature of the heavily obscured 
power source, optically thin tracers in the X-ray, IR, or ra- 
dio wavelengths are necessary. Here, we examine the proper- 
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ties of AzTEC sources identified in the GOODS-South and 
GOODS-North fields using several well-established diagnos- 
tic tests utilizing these optically thin tracers. 

6.1 Spitzer IRAC color-color diagram 

The first set of diagnostic tests to examine are the IRAC 
color-color diagrams that are commonly used for identify- 
ing heavily o bscured AGNs bas ed on the color combinati ons 
proposed by lLacv et al.l l|2004h and IStern et al, I (|2005h as 
shown in Fig. [12] and Fig. 1131 respectively. The majority of 
robustly identified AzTEC sources in both GOODS fields 
fall within the regions previously identified with power-law 
AGNs in both plots. As argued in detail below, these re- 
sults are not a direct consequence of using a Spitzer IRAC 
color selection for identifying AzTEC counterpart sources. 
Interpreting these results should also require some care as a 
young starburst population at z > 1 has a characteristic red 
SED in this part of th e spectrum and should fall within the 
same color region (see lYun et al.|[2008f ). 

The sample size of AzTEC GOODS sources plotted in 
Figs. I12l fc ll3l is nearly 3 times larger than the sample previ- 
ously analyzed by Yun et al., and these new plots show more 
clearly that these AzTEC sources cluster densely around 
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Figure 12. A Ss.Oh/^.S/j versus S5. 8^/53. 6m IRAC color diag- 
nostic diagram for heavily obscured st arbursts and AGN s based 
on the color combinations proposed by lLacv et a l. (2004). Areas 
occupied by dusty young starbursts as noted bv lYun et al. (2008, 
see their Fig. 1) are outlined by a dot-dashed line while the area 
previously identified with power-law spectrum AGNs by Lacy et 
al. is outlined using a dotted line. The IRAC counterparts iden- 
tified with the AzTEC sources in GOODS-South and GOODS- 
North with empty (undetected in X-ray) and dotted circles (de- 
tected in X-ray) cluster around the theoretical dusty starburst 
SED tracks with different amoun ts of dust extincti on (Ay = 
50 fc 200, lEfstathiou et al.l |2000|) . as discussed by lYun et al.l 
(2008). The thick solid line represents the theoretical track ex- 
pected of purely power-law IR AGNs, w hile asterisks are power- 
law A GNs reported by Lacy et al. and iMartinez-Sansigre et all 
(2008). Typical uncertainties in the colors are shown by the cross 
in the upper left corner. 
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Figure 13. A [3.6u] - [4.5/i] versus [5.8 p] - [8.0p] IRAC color 
diagnostic diagram for heavily obscured st arbursts and AGN s 
based on the color co mbinations propos ed bv lStern et al ] {20051 ). 
adopted from Fig. 2 bv lYun et al.l j2008h . All symbols and models 
shown are identical to those in Fig. 1121 



the theoretical col or tracks of 20-80 Myr old dusty star- 
bursts at z > 2 bv lEfstathiou et alJ (|2000h . The dispersion 
in the model tracks and the observed color are larger in 
Fig. 1131 but the AzTEC sources again mostly occupy the re- 
gion spanned by the starburst model tracks, rather than the 
area surrounding the IR power-law track. It is particularly 
noteworthy that AzTEC sources identified with a Chandra 
X-ray source (dotted circles) occupy largely the same area 
as those without X-ray detection (empty circles), and only a 
small fraction of sources (both with and without X-ray de- 
tection) overlap with the IR power-law AGN (asterisks; see 
I Yun et all 120081 ). Conversely, many of the AzTEC sources 
appearing among the power-law IR sources are undetected 
in the 2 fc 4 Ms Chandra surveys ^Alexander et al", I l2003l ; 
iLuo et al.ll200"sl ; I Johnson et ai1l201ll ) . There is a weak trend 
of an increasing fraction of X-ray detected sources with red- 
der color. We can conclude from these diagnostic plots that 
nearly all of the AzTEC sources identified in the GOODS 
fields have IRAC SEDs consistent with that of a young star- 
burst, while a small fraction (< 20%) show IR colors of a 
power-law AGN. 

One thing to clarify is that the use of IRAC colors 
as a method to identify the AzTEC counterparts does not 
lead directly to these observed trends. The adopted color 
selection, [3.6{i] — [4.5/i] ^ 0.0, imposes no restriction on 
the 5.8 U-m and 8.0 urn photometry. This color selection 
is also only one of three independent criteria we examine 
jointly, and the radio data contribute overwhelmingly to the 
secure counterpart identification. In fact, none of the ro- 
bustly identified sources in Table Q] are based on the IRAC 
color selection alone (see § 12.11) . This color selection is not 
used to reject any counterpart candidates either, and in- 
deed several robust counterparts shown in Fig. [T3] have a 
blue IRAC color ([3.6ii] - [4.5/i] < 0.0). In the context of 
the starburst SED model tracks shown, the color selection 
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of [3.6/i] — [4.5/i] ^ 0.0 effectively imposes a redshift bias 
against sources at z < 1 such as AzTEC/GS20. However, 
this bias is more than compensated by the radio and MIPS 
identification methods that systematically favor low-z can- 
didate counterparts. 



6.2 Optical-IRAC-MIPS colors 

Another widely used AGN diagnostic d iagram is the Spitze r 
mid-IR color plot first introduced by llvison et al l (2004). 
They noted that starburst and AGN color tracks as a func- 
tion of redshift are well-separated in the plot of flux ratios 
Su/Sa.o vs. SS.0/S4.5 due to contributions by a power-law 
AGN and PAH emission in the 24 [im band. By analyz- 
ing Spitzer Inf rared Spectrograph (IRS) data on 24 urn se- 
lected galaxies. iPope et al.l (|2008bl ) found that the main dis- 
criminatory information resides in the Sg.o/Si.s flux ratio, 
or the spectral slope in the rest-frame near-IR band. This 
color selecti on is similar to th e ear lier IRAC color an alysis 
proposed bv lLacv et~ai] |2004h and IStern et all (|2005l ). and 
it is again subject to the same confusion with young star- 
burst systems as noted bv lYun et al.l (120081 ) and others. The 
new AGN diagnost ic condition of Ss.o/Si.s ^ 2 proposed by 
IPope et all (|2008bl 'l corresponds to ^(Ss.o/St.s) ^ +0.3 in 
Fig. 1121 and objects satisfying this criteria should also in- 
clude dusty young starbursts at z ^ 0.5 and z > 3 as well 

as power-law IR AGNs. 

Our version of the lPope et all l|2008bl ) diagnostic test is 
shown in Fig. 1141 which shows the plot of the MIPS 24 um 
to optical i-band flux ratio versus Ss.om/^.Sm IRAC band 
flux ratio for AzTEC sources in both GOODS fields. Among 
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Figure 14. A diagnostic color-color diagram using MIPS 24 |xm 
to optical i-band flux ratio versus Ss.0;i/S4.5 M IRAC band flux 
ratio. Symbols are identical to those in Fig. 1121 The long-dashed 
line at S8.0/S4.5 = 2 is the dividing lin e for AGNs and star- 
bursts as proposed by I Pope et al. (2008b). The dotted line near 
lg(5"24/5i)=2.6 is the equivalent division line for the 2 ~ 2 "dust- 
obscu r ed galaxies" with F u {2Afj,rn) / F u (R) > 1000 - see lDev et al.l 
j200Sft : lFiore et al.l <200St) . 



the 57 sources plotted, only 16 (28%) have the Sg.o/Si.s flux 
ratio consistent with hosting an energetic AGN (right of the 
long-dashed line). Again, X-ray detection (dotted circles) 
appears to have little bearing on whether an object falls on 
the starburst ("SB") side or the "AGN" side. The source 
with the highest ratio Ss.cW&.Sp = 4.6 is AzTEC/GS20, 
which is a z — 0.0369 galaxy whose S's.Om/^.s^ ratio arises 
from the bright PAH line emission in the 8 [im band, rather 
than due to a power-law AGN - also see Fig. Q2] 

The choice of the flux ratio between the MIPS 24 urn 
and the optical i— band for the vertical axis in Fig. [14] is 
motivated by the claim of a new class of faint MIPS 24 um 
sources that were missed by earlier optical s tudies. These so- 
called "dust-obsc ured galaxies" ("DOGs"; iDev et alj|2008l : 
iFiore et al.l 12008? ) represent a population of infrared bright 
galaxies that are extremely faint in the optical bands, char- 
acterized by F„(24:fim)/F V (R) > 1000. These z ~ 2 galaxies 
have similar projected density as SMGs and may account for 
as much as ~l/4 of the IR luminosit y density at this redshift 
jDev et al.ll2008l ; |Pope et alj|2008bh . Evidence for both star 
formation and AGN activity has been reported for these ob- 
jects. Based on th e ir sta cking analysis of the X-ray hardness 
ratio, IFiore et alj l|2008h proposed that as many as 80% of 
these dust-obscured galaxies host a Compton-thick AGN. 

Adopting a mean color o f r — i ~ +1 .0 for 
the radio-selected SMGs (e.g., llvison et al.l [2002), the 
u F l/ (24:fim)/F l/ (R) > 1000" definition for DOGs translates 
to l 'lg(S 2 ^/Si) > 2.6" in Fig. M Among the 12 X-ray 
detected secure AzTEC counterpart sources with sufficient 
optical and Spitzer data to be included in this analysis, 
nearly equal numbers of SMGs fall on either side of this 
division line. One significant difference is that 9 X-ray de- 
tected AzTEC sources fall on the optically faint side (above 
the dotted line), while only 3 X-ray detected AzTEC sources 
(including the low-z source AzTEC/GS20) are found on the 
optically bright side. This trend is consistent with the sug- 
gestion by Fiore et al. that many of these dust-obscured 
galaxies host a Compton-thick AGN. However, there are just 
as many AzTEC sources undetected in the X-ray above the 
division line. Given that the observed X-ray emission can be 
largely account ed for by the starburst activity in many cases 
(see below and I Alexander et alj|200"3 ). the significance of a 
higher frequency of X-ray detection among these optically 
faint SMGs is not entirely clear. 

6.3 X-ray 

The 2 fc 4 Ms Chandra X-ray surveys of the two GOODS 
fields ^Alexander et all 120031 iLuo et all 120081 : IXue et all 
2011; Johnson et al.ll201ll ) are some of the deepest X-ray 
data available and thus offers the best opportunity to deter- 
mine X-ray properties of all types of extragalactic sources. 
A cross-correlation of the AzTEC and Chandra X-ray cat- 
alogue has shown that 16 (out of 48) and 8 (out of 40) 
AzTEC sources in the GOODS-South and North fields have 
an X- ray source within 6 " of the AzTEC centroid posi- 
tions (j Johnson et al.ll201lT ). Given the low density of X-ray 
sources, these coincidence are highly statistically significant. 
On the other hand, only a subset of these X-ray sources 
are robust counterparts, and another physical link such a s 
the clustering of massive galaxies (e.g.. lAlmaini et alj|2003h 
must play a role. 
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The derived 2-10 keV X-ray luminosity of these Chandra 
sources favor the AGN-origin for the observed X-ray emis- 
sion. Examining the X-ray proper ties of faint r adio s ources 
in the Hubble Deep Field (North) . lBauer et alj (|2002h found 
that the linear correlation between X-ray luminosity and 
1.4 GHz radio luminosity density of late type galaxies ex- 
tends to luminous X-ray detected emission line galaxies at 
intermediate redshift, suggesting both the X-ray and ra- 
dio processe s are a ssociated with star formation activities. 
IPersic et al.l |2004l ) have shown that the integrated emis- 
sion from high-mass X-ray binaries (HMXBs) can offer a 
natural explanation for the observed correlation, and given 
their short lifetime the measured X-ray luminosity can of- 
fer an instantaneous snapshot of the ongoing star formation 
rate. Since HMXBs also display a characteristic hard X-ray 
spectrum, the hardness ratio of the observed X-ray emis- 
sion does not pro vide a unique probe of AGN activity (e.g., 
iFiore et alj|2008l . - see §E2J. On the other hand, X-ray lu- 
minosity of these Chandra sources associated with AzTEC 
detection ranges between L x (2.0- 10 keV) = 10 42 and 10 43 
erg s -1 . When c onver ted to a SFR using the relation given 
by IPersic etafl (|2004l ), their inferred SFR ranges between 
10 3 and 1O 4 M0 yr _1 , exceeding the SFR derived from their 
UV and IR properties (see Fig. [SJ. A fainter X-ray source 
with 2-10 keV luminosity of 10 42 ergs s _1 can be either a low 
luminosity AGN or an SMG with a SFR = 1O 3 M yr -1 , 
but any source with a higher X-ray lu minosity would requ ire 
a significant AGN contribution (see Ijohnson et al.|[201ll) . 

One intriguing trend found is that the fraction of robust 
AzTEC counterparts that are also CTiandra-detected X-ray 
sources is higher for the brighter AzTEC sources. The 30% 
(7 and 8 out of 25) of the brightest AzTEC sources in the 
GOODS-South and GOODS-North fields are detected indi- 
vidually as a Chandra X-ray source, and an increasing AGN 
activity may be associated with the most luminous AzTEC 
sources. In the same vein, the number of candidate Chandra 
counterpart in the GOODS-South field did not change from 
the 2 Ms catalog to 4 Ms catalog, and the greater depth 
of the X-ray data had curiously little impact. The higher 
frequency of X-ray counterpart in the GOODS-South field 
(16/48 vs. 8/40) may reflect the cosmic variance in these two 
relatively small size fields. A more detailed discussion of the 
X-ray properties of AzTEC sources is presented elsewhere 
l| Johnson et al.ll201ll ). 

The high detection rate for the AzTEC sources in 
the X-ray bodes well for coeval mass growth scenarios for 
the stellar component and the central supermassive black 
hole (SMBH) designed to explain the apparent correla- 
tion between the central blackhole mass and stellar ve- 
locity dispersion ("M-g relation" - Magorrian et alj Il998l ; 
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pie, thr ough detailed numerical modeling, iNaravanan et all 
(|2010bl lar) have shown that a rapid build-up of stellar mass 
and the growth of the central supermassive black hole can 
be achieved through a merger-driven starburst, and can re- 
produce the observed properties of SMGs and dust-obscured 
QSOs. Winds driven by the starburst and the AGN activity 
can effectively disrupt the central concentration of gas and 
dust, driving the evolu tion of such objects fr om an SMG 
phase to a QSO phase (|Naravanan et al1l2008l ). 

A natural consequence of such a scenario is that a mas- 
sive stellar galaxy with a maturing young stellar popula- 



tion would emerge unobscured as the central AGN begins 
to dominate the overall energetics. As the feedback process 
starts to clear out the obscuring dust and gas, the central 
AGN would also become more detectable in the X-ray, UV, 
and optical bands, marking the beginning of the classical 
QSO phase. However, the X-ray detected AzTEC sources 
in the optical and near-IR bands span the entire observed 
range of brightness, indistinguishable from the X-ray un- 
detected sources in Fig. [5] The i — K colors of the X-ray 
detected sources are also indistinguishable from the others, 
suggesting that either (1) the X-ray detection does not signal 
the emergence of the central AGN as the dominant energy 
source or (2) additional complexity is required in the SMG- 
QSO evolution model. 

By analyzing the properties of a sampl e of z ~ 2 SMGs 
exhibi ting broad Ha and H/3 emission lines. [Alexander et al.l 
(2008) have estimated their black hole mass to be > 3 
times smaller than those found in comparable mass nor- 
mal galaxies in the local universe, and > 10 times smaller 
than those predicted for z w 2 luminous quasars and 
radio galaxies. Based on this evidence, they argued that 
the growth of the black hole lags that of the host galaxy 
in SMGs. We find only marginal evidence for AGN con- 
tribution to the near-IR (Figs. [T5] & I13p and mid-IR 
(Fig. I14[) SEDs for the robust AzTEC sources, even among 
those detected in the X-ray. The spectral decomposition 
and the analysis of the emission, absorption, and con- 
tinuum fea t ures in the Spitzer IRS spectra of SMG s by 
Pope et alj (1200881). iMenendez-Delmestre et all \200$ ), and 
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the luminosity in the large majority of cases, even when the 
sample is selected to have AGN-like colors (jCoppin et al.l 
l2010bl ). A consistent trend emerging from these multiwave- 
length data analyses is that dust-obscured starburst activity 
can account for most of the luminosity in SMGs, with little 
or only a minor contribution from AGNs. 



7 SUMMARY AND DISCUSSION 

Taking advantage of some of the deepest imaging and pho- 
tometry data and extensive spectroscopic information in the 
GOODS-South field, we searched for counterparts to the 
48 AzTE C sources fo u nd in the deep 1.1mm wavelength 
survey by IScott et al.l (|2010l ). using a P-statistic analysis 
involving VLA 1.4 GHz, Spitzer/MIPS 24/j.m, and IRAC 
catalogs, combined with cross identification with LABOCA 
870 /J,m sources. Robust (P ^ 0.05) and tentative (0.05 < 
P ^ 0.20) counterpart candidates are found for 27 and 
14 AzTEC GOODS-South sources, respectively. Five of the 
sources (10%) have two robust counterparts, supporting the 
idea that these SMGs are strongly clustered and/or con- 
fused. A spectroscopic redshift is available for 12 robust 
counterparts and 12 tentative counterparts while photomet- 
ric redshifts based on rest-frame UV-to-optical and radio- 
millimetric SED analysis are available for the remainder. 
Stellar mass (M*) and SFRuv are derived by modeling 
the observed optical and Spitzer IRAC photometry while 
SFRir is derived by analyzing the IR, (sub)millimetre, 
and radio photometry using theoretical templates. To im- 
prove the statistics of the subsequent analysis, we applied 
the same counterpart identification and SED analysis to the 
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AzTEC 1.1mm sources identified in the GOODS-North field 
jPerera et al.ll2008l ; iDownes et afl|201l[ ). 

Estimates of the redshift distribution of AzTEC- 
GOODS sources are constructed by combining spectroscopic 
redshifts with UV+optical and radio-millimetric photomet- 
ric redshifts, and these two redshift distributions agree well 
with each other as shown in Figure [4] Our analysis shows 
that 80% of AzTEC sources are at z ^ 2, with a median red- 
shift of z me d ~ 2.6, and there is a significant high-redshift 
tail with 20% of AzTEC sources at z > 3.3. These distri- 
butions are quite diffe rent from the common ly cited SMG 
redshift distribution of lChapman et al.1 l|2005f ). primarily at 
the low redshi ft end. Th e SHA DES survey redshift distribu- 
tion by I Aretxaga et all ([2007) is in better agreement with 
our redshift distribution derived from the AzTEC GOODS 
surveys, and like ours, is missing the low-redshift tail seen in 
IChapman et al.l ()2005h . Complete CO spectroscopic redshift 
surveys using the LMT and ALMA will be able to accurately 
determine the SMG redshift distribution by overcoming the 
large number of systematic biases inherent in all of these 
analyses. 

An examination of the rest-frame UV and optical pho- 
tometry for the securely identified AzTEC sources shows a 
nearly 10 magnitude (a factor of 10 4 in flux density) spread 
in the i— and A"— band photometry and extremely red col- 
ors spanning i — K color between and +6. There are a 
small minority of SMGs that are bright in the rest frame 
UV bands, overlapping with star forming galaxy population 
previously identified in the rest-frame UV searches. On the 
other hand, AzTEC GOODS sources are on average quite 
red and faint, with a median brightness of i = 25.3 and 
K = 22.6, and a large fraction of AzTEC sources are en- 
tirely missed by previous surveys of star forming galaxies. 
Examining the observed scatter in the i — K color, we deduce 
that there is at least a factor of 10 variation in the intrinsic 
rest frame optical luminosity among these SMGs. 

A Hubble diagram of the observed IRAC 3.6 \xm flux 
density shows that these AzTEC-GOODS sources are some 
of the most luminous galaxies in the rest-frame optical bands 
at z ^ 2, offering a good explanation as to why nearly ev- 
ery SMG identified with interferometric observations shows 
a relatively bright IRAC counterpart. Modeling of the ob- 
served rest-frame UV and optical SEDs shows that the stel- 
lar masses are rather large, M* = (1-30) x10 10 Mq, with a 
surprisingly large SFR UV > 100 - 1000 M yr -1 for about 
50% of these galaxies. In comparison, SFRm derived from 
modeling the IR to radio SED covers a relatively tight range 
of 200-2000 Mq yr -1 , independent of the redshift or stellar 
mass. Whether a galaxy has been detected in the X-ray by 
Chandra appears to have no influence on the derived M*, 
SFRuv, SFRm, and SFRuv/SFRih ratio, and the pres- 
ence of an X-ray bright AGN appears to have relatively little 
influence on these quantities. 

These AzTEC-GOODS sources have a specific star for- 
mation rate 10-100 times higher (SSFR « 1-100 Gyr" 1 ) 
than similar stellar mass galaxies at z = 0, and they extend 
the previously ob served rapid rise o f SSFR with redshift 
(SSFR oc (l + z) 5 . [Pamen et alJl^OOSh to z = 2-5. More im- 
portantly, all of the AzTEC-GOODS sources have a SSFR 
above the inverse Hubble time line, indicating that they have 
a current SFR high enough to have built up their entire stel- 
lar mass within the Hubble time at their observed redshift. 



This might be the best quantitative evidence yet that we are 
witnessing these galaxies during their rapid mass build-up 
phase. The flat SSFR as a function of redshift we deduce 
contradicts the model prediction of a tight M*-SFR relation 
based on cosmological h ydrodynamic simulatio ns incorpo- 
rating galactic outflows (|Dave et al.ll2010l . l201ll ). However, 
AzTEC sources with M* > 1O 1O ' 5 M appear to follow this 
model prediction, and one plausible explanation is that a 
different mechanism is operating for the lower mass SMGs, 
leading to a 10 times higher SSFR. Alternatively, they are 
seen at a different phase of rapid mass build-up. Much of the 
cosmic IR background (CIRB) is expected to be associated 
with fainter sources below the confusion limit of our AzTEC 
surveys, and their number counts can offer an important 
constrain to the evolution model for SMGs |Granato et al.l 
12004 iBaugh et al.ll2005l ; H owan- Robinson! I2OO9I ). It still re- 
mains to be shown whether these are young, lower mass 
galaxies seen in their rapid formation epoch or the simmer- 
ing activities in more massive galaxies that have already 
undergone an SMG-like rapid build-up phase in even earlier 
epochs. 

Lastly, we examine the evidence for luminous AGNs in 
these systems using three different diagnostic tests: (1) the 
Spitzer IRAC color-color diagram; (2) optical-IRAC-MIPS 
colors; and (3) X-ray luminosity. We find only marginal 
evidence for AGN contribution to the near-IR (Figs. [12] 
& QUI and mid-IR (Fig. SEDs for the robust AzTEC 
sources, even among those detected in the X-ray. A consis- 
tent trend emerging from this multiwavelength data analysis 
and similar studies by other groups is that dust-obscured 
starburst activity can account for most of the luminosity 
in submm/mm-selected galaxies, with little or only a minor 
contribution from AGNs. 
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APPENDIX A: NOTES ON INDIVIDUAL 
SOURCES 

AzTEC/GSl. There is one clear, robust counterpart 
(GSla in Table [2]), which is a radio source found 4.8" north 
of AzTEC/GSl (Pi.4 = 0.045). This faint radio source has 
an IRAC/MIPS counterpart (P 2 4 M = 0.161), which is also 
a Chandra/ X-ray source. This source has a red IRAC color 
with [3.6/im] — [4.5/im] = +0.37 (P co ior = 0.133), similar to 
the AzT EC sources identified using submillim etre interfer- 
ometry (|Younger et al.ll2007i ; lYun et alJl200Sl ). The 870 pm 
LABOCA source LESS J033211.3-275210 (Ssro u = 9.2+1.2 
mjy) position is only 2.1" away from GSla, and lBiggs et al.l 
(2011) also identify GSla as the robust counterpart. No 
spectroscopic redshift is available for this extremely faint 
optical source (i > 23.5, K > 22.9), and it is likely a 
high redshift (z ^ 3) source with z opt = 2.96 ± 0.45 and 
zmr = 3.56t?jo- 

AzTEC/GS2.1. There are two radio sources (GS2.1a & 
GS2.1b) within the 4.5" search radius, each with a high like- 
lihood of being the AzTEC counterpart. The western- most 
source GS2.1a with a red IRAC color ([3.6/im] — [4.5/im] = 
+0.38) is a robust identification (P 1A = 0.001). The 870 um 
LABOCA source LESS J033219. 0—275219 (S 8T0f i = 9.1+1.2 
mjy) position is only 4.8" away from GS2.1a. The adjacent, 
second faint radio source GS2.1b is an extremely rare ex- 
ample of a faint radio source without any counterpart in 
the IRAC and MIPS images. Both radio sources are for- 
m ally consider e d rob ust counterparts by our analysis and 
bv lBiggs et all (|201lD . There are two additional faint radio 
sources just outside the search radius, making this an ex- 
ceptionally crowded field in the radio band. These two more 
distant radio sources have the same spectroscopic re dshift 
of z = 1.097 (|Le Fevre et alj|2004 iNorris et afll2006l ). and 
their blue IRAC color, [3.6/im] — [4.5/im] = —0.27 suggests 
that they are indeed foreground sour ces. The catalog po si- 
tion for the BLAST 250 um source 59 l|Dunlop et alj|2010h is 
located near the peak of the AzTEC/GS2 contours but be- 
tween the two deconvolved components AzTEC/GS2.1 and 
AzTEC/GS2.2. A third potential counterpart, GS2.1c, iden- 
tified by MIPS 24 um detection is only 3.9"from the AzTEC 
centroid; however it has a blue IRAC color ([3.6/im] — 
[4.5/im] = —0.32) with a spectroscopic redshift of z — 0.644, 
and is likely a foreground object. 

AzTEC/GS2.2. There are no radio sources within 15" of 
AzTEC/GS2.2. The MIPS 24 pm source GS2.2a is a po- 
tential counterpart with Pia^l = 0.212 with a blue IRAC 
color ([3.6/im] — [4.5/i m] = —0.23) and a spectroscopic 
redshift of z = 1.046 (jPopesso et alj|2009h . Therefore, it 
is likely a member of the foreground z — 1.10 cluster 



GCL J0332. 2-2752 {a v = 433 km/s; iDiaz-Sanchez et all 
l2007h whose center is located only ~ 20" to the north- 
east, at a = 03 h 32 m 17.5" and 5 = -27°52'32". This MIPS 
source is blended with a second IRAC source located 3" 
to the southwest, GS2.2b, which has a red IRAC color 
([3.6/im] - [4.5/im] = +0.26; P coior = 0.390). GS2.2b is also 
a BzK galaxy and thus is an actively star forming system. 
The IRAC source GS2.2c is another BzK galaxy with a red 
IRAC color ([3.6/im] — [4.5/im] = +0.10) and is an inter- 
esting alternative counterpart, though P co ior = 0.439. The 
closest 870 pm LABOCA source LESS J033217.6-275230 
(•SWop — 6.3 ± 1.3 mjy) position is 15" northeast of the 
AzTEC centroid, nearly centered on the z = 1.1 cluster. 
The position of AzTEC/GS2.2, however, is very uncertain 
as it is blended with AzTEC/GS2.1, so the counterpart iden- 
tification is highly problematic. 

AzTEC/GS3. The faint IRAC source GS3a (Siaghz = 
40.7 ± 6.5 /iJy; P 1A = 0.045) is also a MIPS 24 pm 
source with a red IRAC color ([3.6/im] — [4.5/im] = +0.37; 
Pcoior = 0.174). The 870 pm LABOCA source LESS 
J033248. 1-275414 (S 870m = 8.8 ± 1.2 mjy| is only 2.4" 
away from GS3a, leading iBiggs et alj (|201ll) to conclude 
this s ource as a secure counterpart as well. iDunlop et al.l 
also identified GS3a as the counterpart to the 250 
pm BLAST source 593 and estimated a redshift z > 2.5 for 
this optically invisible source. Our radio-mm photometric 
redshift of zmr = 3.09jlx"ii supports this high-z hypothesis. 
There is a faint red IRAC source, GS3b, which is is also a 
tentative detection that cannot be ruled out. 

AzTEC/GS4. There is only one red IRAC source (GS4a; 
Pcoior = 0.070) within the 6.5" search radius. It is a faint ra- 
dio emitter (Siaghz = 25.4 ± 6.5 /iJy, Pi a = 0.021). There 
are two other radio sources within 17" from the AzTEC posi- 
tion, but GS4a is the only source falling within the error cir- 
cle of the 870 pm LABOCA source LESS J033249. 2-274246 
(SWop = 8-8 =t 1-2 rnjy). Therefore GS4a is a robust coun- 
terpart for the AzTEC source although it is only a ten- 
tativ e counterpart for the LABOCA source (IBiggs et al.l 
l201lf ). This is another high-z candidate source with z op t = 
3.37 ± 0.25 and zmr = 3.53t?;27. 

AzTEC/GS5. There is a single tentative counterpart 
within the 7.1" search radius from the AzTEC/GS5 posi- 
tion. However, the AzTEC contours are elongated in the 
east-west direction, joining the two VLA radio sources GS5a 
& GS5b. Both radio sources have red IRAC counterparts, 
and both sources may contribute to the AzTEC emission. 
The 870 pm LABOCA source LESS J033150.8-274438 
(<S , 870m — 3.9 ± 1.4 mjy) is located on ly 4.1" away from 
GS5 a (also a Chandra X- ray source), and lBiggs et al.l (2011) 
and IChapin et al.l (|201ll ) also identi fy GS5a as t he se cure 
counterpart to the LABOCA source. ICasev et all i|201ll ) re- 
ported a spectroscopic redshift of z = 1.599 for GS5a based 
on some absorptions features, and this redshift is further 
supported by the 9 h r long integration VLT spectrum by 
ISilverman et all (|2010h . 

AzTEC/GS6. The IRAC/MIPS source GS6b is located 
only 0.4" from the AzTEC/GS6 centroid. However, its blue 
IRAC color ([3.6/im] — [4.5/im] = —0.23) and spectroscopic 
redshift of z — 1.102 (Stern et al., in prep) suggest that 
GS6b is likely a foreground object. At a distance of 12.4", 
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the IRAC/MIPS source GS6a is located outside the 7.5" for- 
mal search radius for a counterpart, but it has very red IRAC 
color ([3.6/im] — [4.5/im] = +0.45) and is a faint radio source 
{S 1A ghz = 31.0±6.3 /ijy; P XA = 0.268). iBiggs et alj <|201lT ) 
identify GS6a as the robust counterpart for the 870 [im 
LABOCA source LESS J033225.7-275228 (SWom = 5.8±1.4 
mjy) located 6.2" away, and this galaxy is likely the primary 
counterpart to the AzTEC source as well. 

AzTEC/GS7. The red IRAC/MIPS source GS7a is the 
only radio source within the 8.7" search radius (Siaghz = 
51.2 ± 6.4 /iJy; Pi. 4 = 0.126). Therefore, it is considered 
a tentative counterpart to AzTEC/GS7. The spectroscopic 
redshift of GS7a, which is also identified as a Chandra X- 
ray source, is z — 2.676, in excellent agreement with its 
radio -mm photometric redshift of zmr — 2.56^ 092- 

The 

870 um LABOCA source LESS J033213.6-2756 02 (S S70u = 
9.1 ± 1.2 mjy) is located only 4.2" away, and Big gs et al.l 
l|201lf) also identify GS7a as a robust counterpart. A sec- 
ond possible counterpart, GS7b, is a red IRAC/MIPS source 
(p24/j = 0.151, Pcoior = 0.168) located on the other side of 
the AzTEC centroid and may contribute to the observed 
1100 \im emission. 

AzTEC/GS8. There is a single robust radio counterpart 
GS8a, 4.4" from the AzTEC centroid position with Pi. 4 = 
0.038. The IRAC/MIPS counterpart GS8a also has a red 
IRAC color ([3.6/im] — [4.5/im] = +0.33) and relatively 
bright MIPS 24 urn emiss i on (S 2 4 U = 620 /ijy; P 2Al , = 
0.203). Both IChapin et all (|201ll ) and iBiggs et all (|201lh 
have identified this z — 2.252 galaxy as the counterpart 
to the 870 nm LABOCA source LESS J033205. 1-274652 
(SWcv = 7.5+1.2 mjy), located only 7.7" away from AzTEC 
position. The second IRAC/MIPS source GS8b, located 2.9" 
from the AzTEC position, is also a plausible MIPS 24 um 
counterpart (p24 M = 0.089) with red IRAC color and is also 
a Chandra X-ray source. 

AzTEC/GS9. The single tentative radio counterpart 
(Si. 4 = 86.8 ± 6.6 /iJy; Pi. 4 = 0.070) is associated with 
a red IRAC/MIPS source GS9a, located 6.2" from the 
AzTEC centroid position. It is also a CTiandra-detected X- 
ray source and should be considered a strong candidate for 
the AzTEC counterpart. Located only 3.5" away from the 
AzTEC centroid, the IRAC/MIPS source GS9b is an in- 
triguing alternate counterpart candidate given its red IRAC 
color ([3.6/im] - [4.5/^m] = +0.29) and MIPS 24 um emis- 
sion. If GS9b is the primary source of 1100 (xm continuum 
emission, then it is likely to be a high redshift system as 
its radio non-detection requires zmr > 3.3. Slightly further 
away is GS9c, also a tentative red IRAC source. No nearby 
source is found in the 870 \im LABOCA catalog, but the 
LABOCA map shows a S/N = 3.0 peak (3.4 ± 1.2 mjy) at 
the position of GS9a. 

AzTEC/GSlO. There is a single robust radio counter- 
part to AzTEC/GSlO, located 5.3" from the AzTEC cen- 
troid (Sia = 89.3 ± 6.4 /iJy; P1.4 = 0.053). Its red 
IRAC/MIPS counterpart GSlOa ([3.6/im] - [4.5/im] = 
+0.14) has a reported spectroscopic redshift of z = 0.0338 
in the GOODS/ESO VIMOS DR1 catalog jPopesso et all 
l2009l ). but the s ame group revised it s redshift to z = 2.035 
using new data (|Balestra et al.ll2010l ). This revised spectro- 
scopic redshift is in excellent agreement with our photomet- 



ric redshift, zmr — 2.03^5) '75- Another red IRAC source 
GSlOb ([3.6/xm] - [4.5/tm] = +0.04; P coior = 0.077), located 
only 3.2" away from the AzTEC position, is not detected in 
the radio or by MIPS. The 870 um LABOCA source LESS 
J033207.6-275123 (S 870 » = 7.6 ± 1.3 mjy) is located only 
4.3" away from GSlOa, a nd it is also i dentifi ed as a robust 
LABOCA counterpart bv lBiggs etafl (|201 lh . 

AzTEC/GSll. There is a single tentative VLA radio coun- 
terpart for this source located 6.5" from the AzTEC centroid 
(S1.4 = 46.0 ± 6.4 /ijy; Pi. 4 = 0.081), and it is also a Chan- 
dra X-ray source. Its IRAC/MIPS counterpart GSlla has a 
slightly blue color ([3.6/im] - [4.5/im] = -0.02), but the 
VLT ISAAK if-band image l|Retzlaff et al.l |2010| ) shows 
that this IRAC source is a blend of an optically bright 
(i — 21.7) z = 0.246 foreground source and an optically 
faint source second source, which is the radio source. The 
250 ^m BLAST source 109, located ~30" southeast of the 
AzTEC centroid, is undetected at longer wavelength bands 
and is identified wi th a z — 0.124 foreground disk galaxy 
l|Dunlop et al.l [2010+ Our photometric redshift for GSlla, 
zmr = 2.501°'^, is completely inconsistent with this identi- 
fication, and the proposed BLAST counterpart is unlikely to 
be related to the AzTEC source. No nearby source is found 
in the 870 um LABOCA catalog, but the LABOCA map 
shows a S/N ~ 3 peak (3.5 ± 1.2 mjy) at the position of 
GSlla. 

AzTEC/GS12. The most likely counterpart candidate is a 
faint radio, red IRAC source GS12a located 4.0" away from 
the AzTEC position with [3.6/im] - [4.5/im] = +0.10. This 
z = 4.762 galaxy was also identified as the counterpart to the 
LABOCA 870 um survey s ource LESS J033229 3- 275619 
(SWop = 5- 1 ± 1A mJ y) by ICoppin et all (120091 ') based on 
its proximity to the LABOCA position and the presence of a 
~ 3a radio source. R edshifted CO e mission at zco = 4.755 
has been reported bv lCoppin et al.l l|2010al ). lending further 
support for the counterpart identification. 

AzTEC/GS13. There is a high concentration of 
IRAC/MIPS sources with spectroscopic redshifts between 
1.0 and 1.6 in the region surrounding AzTEC/GS13. 
The most likely counterpart for AzTEC/GS13 is a red 
IRAC/MIPS source GS13a ([3.6/im] - [4.5/im] = +0.24), lo- 
cated only 2.1" away from the AzTEC centroid. This source 
is also a faint radio source with no known spectroscopic red- 
shift. A second IRAC/MI PS sources GS13b is a Chandra 
X-ray source at z = 1.039 (jMignoli et al.|[2005l ). but its sta- 
tistical likelihood of being the AzTEC counterpart is lower 
(see Table [T}. A third candidate counterpart, GS13c, is a 
faint radio source, but it has a blue IRAC color ([3.6/im] — 
[4.5/tm] = -0.25) and is therefore likely a foreground source. 
The 250 um BLAST source 193 is located ~25" south of the 
AzTEC position. Although the low density of the AzTEC 
and BLAST sources make the chance coincidence of these 
two sources even at such a substantial separation small, 
few plausible candidates are found within the BLAST posi- 
tion error circle. No nearby source is found in the 870 [im 
LABOCA catalog, but the LABOCA map shows a S/N ~ 3 
peak (3.6 ± 1.2 mjy) near the AzTEC centroid position. 

AzTEC/GS14. There are no radio sources within the 9.0" 
search radius. The only tentative counterpart is a faint, red 
IRAC source GS14a ([3.6/im] - [4.5/im] = +0.04; P color = 
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0.083) located only 3.1" from AzTEC/GS14. Although no 
nearby source is found in the 870 urn LABOCA catalog, the 
LABOCA map shows a S/N ~ 3 peak (3.3±1.2 mjy) nearly 
coincident with this position. The previously unpublished 
redshift of z = 3.640 (Keck LRIS spectrum, H. Spinrad, 
priv. comm.) for GS14a is consistent with its photoz z opt = 
3.50 ±0.30 and the non-detection in the radio and the MIPS 
24 um bands, similar to the z = 4.762 AzTEC/LABOCA 
source GS12a. 

AzTEC/GS15. There are several faint, red IRAC/MIPS 
sources within the 9.0" search radius, although they are indi- 
vidually not particularly compelling. The faint radio source 
GS15a, although 12. 7" from the AzTEC centroid, is also a 
Chandra X-ray source and is located 1.4"from the centroid 
of the LABOCA 870 urn source LESS J033151. 5-274552 
(SWom = 4.2 ± 1.4 mjy), and G S15a is the secure LABOCA 
counterpart l|Biggs et al.ll2011T ). 

AzTEC/GS16. There are two faint radio sources within 
the counterpart search radius of 10.5". The red IRAC/MIPS 
source GS16a ([3.6/im] — [4.5/im] = +0.53) is 6.1" away from 
the AzTEC position, and it is a tentative radio counterpart 
with Pi. 4 = 0.072. GS16a is also an X-ray source detected by 
Chandra and has a spectroscopic redshift of z = 1.719 (Sil- 
verman et al., in prep.). The second radio source GS16b is 
located slightly further away, 7.9". Although it is a brighter 
MIPS 24 urn source, it has a blue IRAC color ([3.6/im] - 
[4.5/im] = —0.29) an d a spectroscopic redshift of z = 1.017 
jMignoli et ail [2005). suggesting it is a foreground source. 
No nearby source is found in the 870 um LABOCA cata- 
log, but the LABOCA map shows a S/JV ~ 2 peak near the 
position of G16a. 

AzTEC/GS17. Two plausible IRAC/MIPS sources are 
found within the 10.5" search radius. The IRAC/MIPS 
source GS17a is a faint radio source with Pi. 4 = 0.007, 
making it formally a robust identification. However, it has 
a blue IRAC color with [3.6/im] - [4.5/im] = -0.27 and 
is thus likely a foreground object(z opt = 1.01 ± 0.10). The 
IRAC/MIPS source GS17b is another robust identification 
based on the MIPS detection (P 2 4 M = 0.026) and red IRAC 
color ([3.6/im] - [4.5/im] = +0.19, P co ior = 0.021), and this 
Chandra dete cted X-ray source has a photometric redshift 
of 2 = 2.66 (|Silverman et alj|20ld ). Though slightly fur- 
ther from the AzTEC centroid, GS17c is an interesting al- 
ternative possibility: this red IRAC/MIPS source ([3.6/im] 
— [4.5/im] = +0.36) has optical colors satisfying the BzK 
criteria for star forming galaxies at high redshifts. There 
is also a close pair of radio sources, GS17d and GS17e, lo- 
cated 12"north of the AzTEC position; though they are out- 
side the nominal search radius, both have red IRAC colours 
and remain plausible counterparts to the AzTEC source. No 
nearby source is found in the 870 um LABOCA catalog, but 
the LABOCA map shows a S/N ~ 3 emission peak between 
GS17b and the two radio sources GS17d & GS17e. 

AzTEC/GS18. There are two faint radio sources within 
the 9.3" search radius of AzTEC/GS18. The closest 
IRAC/MIPS source GS18a is formally a robust counterpart 
owing to its proximity to the AzTEC position and its very 
red IRAC color ([3.6/im] - [4.5/im] = +0.47, P coloT = 0.113). 
Additionally, GS18b is a z — 2.6 88 galaxy which is an X- 
ray source detected by Chandra l|Akivama|[2005l ) and has 



a red IRAC color ([3.6/im] — [4.5/im] = +0.31), making it 
a plausible candidate for the counterpart to AzTEC/GS18. 
A third radio- faint IRAC/MIPS source GS18c is a tenta- 
tive candidate based on its MIPS 24 um emission, but it 
has very blue colour ([3.6/im] — [4.5/im] = —0.39) and is 
likely a foreground object. The 870 urn LABOCA source 
LESS J033243.6-274644 (S 8 70 M = 6.4 ± 1.3 mjy) is lo- 
cated ~8.4" south of the AzTEC position, pu tting GS18a 
betwe en the AzTEC and LABOCA centroids. iBiggs et ail 
|201ll) identified GS18a as a tentative (and only) counter- 
part to the LABOCA source. This extremely faint optical 
source (i — 28.1, K > 24.5) is another high-z candidate with 
z MR = 3.001?;^. 

AzTEC/GS19. A pair of radio sources, GS19a & GS19b 
(with Pi. 4 = 0.102 and Pi. 4 = 0.143, respectively), are 
strongly favored as the counterparts to AzTEC/GS19 when 
all statistical measures are taken into account. They form 
a blended source in the MIPS 24 um band, and both 
have red IRAC colors. It is possible that these two sources 
are physically related and both contribute to the AzTEC 
emission. No spectroscopic redshift is available for either 
sources while their photometric redshifts are quite similar 
(Zopt = 1.83 ± 0.35 & 2.08 ± 0.25). At z = 2.0, their pro- 
jected separation of 3.2" corresponds to 27 kpc. There is 
a third IRAC/MIPS counterpart, GS19c, that is a tenta- 
tively counterpart with a lower likelihood (P24/J = 0.181). 
No nearby source is found in the 870 um LABOCA catalog, 
but the LABOCA map shows a S/N ~ 2.4 emission peak 
near GS19a & GS19b. 

AzTEC/GS20. The AzTEC contours are well centered and 
follow the light profile of the z = 0.0369 galaxy GS20a, 
which is also a bright and well resolved radio, MIPS 24 & 
70 um, and X-ray source. This source is quite blue ([3.6/im] 
— [4.5/im] = —0.43) and is obviously a bright foreground 
galaxy. It is the brightest BLAST 250 urn source within the 
GOODS-South survey field proper, and lDunlop et all (|2010h 
argue that this foreground galaxy is the correct counter- 
part for the BLAST source. It is possible that the 1.1 mm 
emission originates from the cold dust associated with this 
spiral galaxy (as seen in the HST images), but it is diffi- 
cult to reproduce the entire measured spectral energy dis- 
tribution from A = 1 um to 20 cm using a reasonable set 
of assumptions on dust temperature, IR luminosity, and 
star formation rate for this low redshift galaxy. Instead, the 
AzTEC emission may originate from a background object, 
possibly lensed by this foreground galaxy, similar to AzTEC 
J100008.05+022612.2 imaged at high angular resolution by 
lYounger et all \200H \. 

AzTEC/GS21. The single tentative radio counterpart 
(GS21a) has [3.6/im] — [4.5/im] = +0-24 and a spectroscopic 
redshift of z = 1.910 jvanzella et alj 120081 ). The position 
centroid of the BLAST 250 um source 861 is displaced by 
~7" to the northeast, but GS21a is still within the error- 
circle of the BLAST source. No nearby source is found in the 
870 vim LABOCA catalog, but the LABOCA map shows a 
S/N ~ 1.8 emission peak near GS21a. A second tentative 
counterpart GS21b is a red IRAC/MIPS galaxy located just 
2.5"away from the AzTEC centroid. 

AzTEC/GS22. The faint radio source GS22a, located 7.8" 
away from the AzTEC centroid, is the most likely coun- 
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terpart (Si. 4 = 34.6 ± 6.5 ^y). This galaxy is also a red 
IRAC/MIPS source ([3.6/Ltm] - [4.5/xm l = +0.30) and ha s 
a spectroscopic redshift of 2 = 1.794 l)Wuvts et al.l F2009) . 
The 250 pm BLAST source 552 position centr o id is ~10" 
east of the AzTEC position, and IDunlop et all (|2010t ) also 
identified the radio source GS22a (located ~7" away from 
the BLAST position) as the likely counterpart. A second 
red IRAC/MIPS candidate counterpart, GS22b, is closer to 
the AzTEC centroid and is an interesting alternative. No 
nearby source is found in the 870 pm LABOCA catalog, 
but the LABOCA map shows a S/N ~ 3.1 emission peak 
near this red IRAC/MIPS source GS22b. 

AzTEC/GS23. The faint radio source GS23a (Si. 4 = 
23.4 ± 6.5 A*Jy; -Pi. 4 = 0.030) is a robust counterpart, and 
its red IRAC color ([3.6/im] - [4.5/^m] = +0.45) adds to the 
high likelihood of being the correct counterpart. A second 
tentative radio counterpart, GS23b is only slightly further 
from the AzTEC centroid, and it is a red IRAC /MIPS source 
with a spectroscopic redshift of 2 = 2.277 jChapin et al.l 
120111 ). Both GS23a & GS23b are within the beam area 
of the 870 urn LABOCA source LESS J033221. 3-275623 
(S87 0u = 4.7 ± 1. 4 mjy ) an d the 250 urn BLA ST source 
158. lDunlop et all (|2010l ) and lBiggs et~al1 (|201ll ) have iden- 
tified GS23b as the counterpart consisting of "an extremely 
complex faint system" at optical wavelengths. 

AzTEC / GS24. There is a single tentative faint radio coun- 
terpart within the 12.2" counterpart search radius, GS24a. 
This red IRAC/MIPS source is the most likely counterpart 
to AzTEC/GS24. The 250 urn BLAST source 104 straddles 
GS24a and a z = 2.578 type 2 QSO J033235. 78-274916. 82 
(Rigo poulou et al.|[2009l . also detected at MIPS 70 pm), and 
both sources l ikely contribute to the 250 pm continuum. 
IDunlop et all (|2010l ) instead identified the 2 = 0.547 radio 
source located at the edge of the AzTEC and BLAST beam 
based on the radio P— statistic. No nearby source is found 
in the 870 pm LABOCA catalog, but the LABOCA map 
shows a S/N ~ 3 emission peak near the 2 = 2.578 type 2 
IR QSO J033235. 78-274916.82 and a secondary S/N ~ 2 
emission peak near GS24a. The AzTEC contours are elon- 
gated in the north-south direction, and this may be another 
example of a blended source. 

AzTEC/GS25. The red IRAC/MIPS source associated 
with radio emission GS25a is located only 6.8" away from 
the AzTEC centroid position. This galaxy, also detected 
in the X-ra y by Chandra, has a spectroscopic redshift of 
2 = 2.292 jPopesso et al.l |2009| V The 870 pm LABOCA 
source LESS J033246.7-2751 20 (S 870 u = 5.9 ± 1.3 mjy) 
is well-centered on GS25a, and lBiggs et alj (|201ll ) also iden- 
tify the same galaxy as the robust LABOCA counterpart. 

AzTEC/GS26. AzTEC/GS26 has no radio counterpart 
within the 12.2" search radius of the AzTEC position, but 
there are five red IRAC sources. Of these, GS26a ([3.6^im] 
- [4.5/xm] = +0.22; P coior = 0.237) located 5.5" away is the 
most probable counterpart to this AzTEC source. GS26b 
and GS26c are other red IRAC sources located slightly 
further away (~ 7"). No nearby source is found in the 
870 pm LABOCA catalog, but the LABOCA map shows 
a S/N ~ 2.9 emission peak near the 2 = 2.331 MIPS 
galaxy J033216.3-274343.4 with a red IRAC color ([3.6/xm] 



- [4.5/wn] = +0.35), located ~ 7.5" away from the AzTEC 
position. 

AzTEC/GS27. The red IRAC/MIPS source GS27a 
([3.6^im] — [4.5/xm] = +0.38) is located 13.0" away from 
the AzTEC centroid, just within the counterpart search ra- 
dius of 13.0", and is associated with weak radio emission 
(S1.4 = 23.6 ± 6.5 /uJy) from a 2 = 2.577 l|Popesso et all 
2009) galaxy. However, formally this is not a secure iden- 
tification (Pi. 4 > 0.20) owing to its large separation from 
the AzTEC centroid. Still, this is the most plausible coun- 
terpart within the proximity of AzTEC/GS27. No nearby 
source is found in the 870 pm LABOCA catalog, but the 
LABOCA map shows a S/N ~ 3.1 emission peak near this 
red IRAC/MIPS source GS27a. 

AzTEC/GS28. AzTEC/GS28 has no radio source within 
the 13.0" counterpart search radius. The red IRAC sources 
GS28a ([3.6Aim] - [4.5/xm] = +0.55; P coior = 0.073) and 
GS28b ([3.6/mi] - [4.5/j.m] = +0.07; P color = 0.136) are the 
most likely counterparts. No nearby source is found in the 
870 pm LABOCA catalog, but the LABOCA map shows an 
extended source with a S/N ~ 3.0 emission peak near the 
AzTEC source position, and this may be another blended 
source. 

AzTEC/GS29. There are no compelling counterparts to 
AzTEC/GS29. The closest IRAC/MIPS source GS29a (3.8" 
away) is very blue ([3.6/xm] — [4.5/im] = —0.43) and has 
a spectroscopic redshift of 2 = 0.577, so this is likely 
a foreground source. The next closest source GS29b is a 
2 = 2.340 galaxy located 4.6" away, but it is also blue 
([3.6/^m] — [4.5/im] = —0.10). There are two radio sources 
located just outside the 13.0"search radius, but they both 
have relatively high probabilities of being false associations. 
No nearby source is found in the 870 pm LABOCA catalog, 
but the LABOCA map shows an extended source with two 
S/N ~2-3 emission peaks near the red IRAC/MIPS sources 
J033158.65-274516.3 and J033200. 13-274453.2. 

AzTEC/GS30. The brightest MIPS 24 pm source within 
the counterpart search radius of 13.5" (S24H = 459 ± 6/iJy; 
P24/J = 0.44) is also a faint radio source (Si. 4 = 37.2 ± 
6.2/iJy; P1.4 = 0.082) with a red IRAC color ([3.6/zm] - 
[4.5/im] = +0.19). Therefore, GS30a is a tentative but the 
most plausible counterpart for AzTEC/GS30. The spectro- 
scopic redshift of this galaxy is yet unknown. A second ra- 
dio source (GS30b) and a radio-faint IRAC source (GS30c) 
are also tentative counterparts to AzTEC/GS30, both with 
red IRAC colors. No nearby source is found in the 870 pm 
LABOCA catalog, but the LABOCA map shows an ex- 
tended source with a S/N = 2.9 centered on the AzTEC 
centroid position. 

AzTEC/GS31. Both of the two bright MIPS 24 pm 
sources (GS31b) within the 13.6" search radius region are 
associated with faint radio emission. The western source 
GS31a however is closer (2.7") to the AzTEC centroid, mak- 
ing it a robust counterpart candidate (Pi. 4 = 0.015). The 
eastern source GS31b is located 7.9" away and is a tentative 
counterpart (Pi. 4 = 0.118). It is possible that both sources 
contribute to the 1.1 mm emission detected by AzTEC. 
They both have blue IRAC colors ([3.6/im] — [4.5/im] = 
—0.25 and —0.36 for GS31a and GS31b, respectively) how- 
ever, and the true counterpart may be a much fainter source 
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located between or behind these so urces. The spectros copic 
redshift of GS31a is z = 1.843 l|Wuvts et al.ll2009l ). No 
nearby source is found in the 870 \im LABOCA catalog, but 
the LABOCA map shows an elongated north-south struc- 
ture with a S/N — 2.8 peak, similar to the AzTEC source 
morphology, and it is likely another blended source. 

AzTEC/GS32. There are three radio counterparts within 
the 13.5" search radius of the AzTEC centroid, and all 
three are also MIPS 24 um sources. The nearest radio 
source GS32a (Si a = 30.3 ± 6.8 pJy, P 1A = 0.162; and 
S 2 a i , = 371.1 ± 11.7 pJy, P 2 a i , = 0.421), located 9.6" away 
from the AzTEC centroid, is the primary candidate by its 
proximity. The source GS32b is similarly bright in radio and 
MIPS 24um bands, but it is slightly further away, 10.5", 
from the AzTEC centroid position. Both GS32a and GS32b 
are quite blue, however; we include GS32c, a faint radio 
galaxy with red IRAC color, in the catalog as an alterna- 
tive, although it is located 13.4" from the AzTEC centroid. 
No nearby source is found in the 870 um LABOCA catalog, 
and no significant emission peak is found at this location in 
the LABOCA map. 

AzTEC/GS33. The faint IRAC source GS33a located 7.4" 
from the AzTEC centroid position is associated with a weak 
radio source (5*1.4 = 28.6 ± 6.2 pJy). It is also a weak MIPS 
24 um source and has a red IRAC color ([3.6/im] — [4.5pm] 
= +0.34; Pcoior = 0.390). Therefore this galaxy is a tentative 
candidate counterpart for AzTEC/GS33. No nearby source 
is found in the 870 um LABOCA catalog, but the LABOCA 
map shows an isolated emission peak with S/N ~ 2.3 cen- 
tered near GS33a. 

AzTEC/GS34. There is a high density of IRAC sources in 
this field, but few are located in the immediate vicinity of 
the AzTEC centroid position. There are no robust or tenta- 
tive counterparts to AzTEC/GS34: P > 0.20 for all sources 
within the 13.5"search radius. No nearby source is found 
in the 870 um LABOCA catalog, but the LABOCA map 
shows an emission peak with a S/N ~ 3.5 centered near the 
z — 1.356 (Silverman et al., in prep) faint radio and MIPS 
source GS34a. 

AzTEC/GS35. The red IRAC source GS35a, located 2.0" 
away from the AzTEC source centroid position, has a robust 
radio counterpart (S1.4 = 41.3 ± 6.7pJy; P\ A = 0.008). Its 
MIPS emission (P 24m = 0.027) and red IRAC color ([3.6pm] 
— [4.5/im] = +0.37; Pcoior = 0.022) make this galaxy a 
robust counterpart to the AzTEC source. Its spectroscopic 
redshift is unknown. The second red IRAC source GS35b lo- 
cated 4.9" away from the AzTEC position (P co ior = 0.191) is 
undetected at radio or MIPS 24 um bands. No nearby source 
is found in the 870 urn LABOCA catalog, and the LABOCA 
map shows an emission peak with a S/N ~ 2.1 centered near 
the z = 0.734 radio-loud QSO J033227.00-274105.0. 

AzTEC/GS36. No radio source is detected within the 
13.5" counterpart candidate search radius, and the galaxy 
GS36a is the only IRAC source with a red IRAC color, 
[3.6pm] - [4.5pm] = +0.68 (P coior = 0.304). So this IRAC 
source is the most promising counterpart candidate, al- 
though it has a high probability of false association. Few 
other potential candidates are present in this field. No 
nearby source is found in the 870 um LABOCA catalog, and 



the LABOCA map does not show any significant emission 
peak within the search radius. 

AzTEC/GS37. The red IRAC/MIPS source GS37a with 
[3.6pm] — [4.5pm] = +0.08 is associated with faint radio 
emission (Sia = 20.5 ± 6.4pJy), though it is not formally 
a tentative counterpart candidate for AzTEC/GS37. There 
is a single tentative red IRAC source 4.3" from the AzTEC 
centroid, which is possibly at high redshift (z > 3) given its 
non-detection at 1.4 GHz. No nearby source is found in the 
870 um LABOCA catalog, but the LABOCA map shows a 
ridge of 2.1-2.7cr peaks running along the similar ridge found 
in the AzTEC map (including GS37a), suggesting this is 
likely another blended source. 

AzTEC/GS38. There is a single tentative radio coun- 
terpart (GS38a; Pi. 4 = 0.116), which is also a Chandra 
X-ray source. It has a very blue IRAC color ([3.6pm] — 
[4.5pm] = —0.55) and a spectroscopic redshift of z = 0.735 
IVanzella et~aH 120081 ). Therefore GS38a is likely a fore- 
ground object. No nearby source is found in the 870 
LABOCA catalog, but the LABOCA map shows a 2.3<r peak 
near the red IRAC/MIPS source GS38b, located 7.4" away 
from the AzTEC centroid (P 24fl = 0.453, P color = 0.386). 

AzTEC/GS39. There is a single tentative radio coun- 
terpart (GS39a; Pi. 4 = 0.083) which is also a red IRAC 
source ([3.6/im] — [4.5pm] = +0.55). The position centroid 
of the 870 um LABOCA source LESS J033154.4-274525 
(Ssro u = 3.8 ± 1.4 m jy) is offset from GS39a by only 
4.9". iBiggs et all (|201 ll ) also identify GS39a as a tentative 
LABOCA counterpart, and GS39a is a strong counterpart 
candidate for AzTEC/GS39. Another tentative counterpart, 
GS39b, is a MIPS source located only 2.8"from the AzTEC 
centroid, but given its blue color, it is likely a foreground 
object. 

AzTEC/GS40. There are no compelling counterparts to 
AzTEC/GS40. We list a single IRAC source within the 
15.0"search radius. GS40a is 10. 2" away, and is slightly blue 
with a high probability of being a false association. No 
nearby source is found in the 870 um LABOCA catalog, and 
the LABOCA map shows only a 2a emission peak located 
~ 10" northwest of the AzTEC position. 

AzTEC/GS41. The radio sources GS41a and GS41b are 
both promising counterparts to AzTEC/GS41 with Pi. 4 = 
0.120 & 0.132. They are also red IRAC sources with [3.6pm] 
— [4.5pm] = +0.23 & +0.30 although they are just far 
enough away from the AzTEC centroid to make them 
tentative counterparts only. Two other red IRAC sources 
GS41c and GS41d are slightly closer to the AzTEC cen- 
troid and are tentative counterparts as well. The centroid 
of the 870 um LABOCA source LESS J033302. 5-275643 
(5*870^ = 1 2.0 ± 1.2 mjy) is l ocated closest to GS41b and 
GS41d, and lBiggs et all (|201lh identify GS41b as the robust 
counterpart to the LABOCA source. 

AzTEC/GS42. There are no radio counterparts within 
the 6.9" search radius of AzTEC/GS42. The only tentative 
counterpart is GS42a, a red IRAC source ([3.6pm] — [4.5pm] 
= +0.04) that is also well-matched to the position of the 
LABOCA source LESS J033314.3-275611 (S 8 7o M = 14.5 ± 
1.2 mjy). The AzT EC contours are ex tended in the north- 
south direction, and I Weiss et al.l (|2009h modeled this source 
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with a second, fainter co mponent (LESS J0 33313. 0—275556, 
SW = 4.3 ± 1.4 mjy). IBiggs et al.1 (|201lh identified adja- 
cent red IRAC/MIPS source J033314.41-275612.0 as the 
robust LABOCA counterpart. 

AzTEC/GS43. The only tentative counterpart is a red 
1RAC source GS43a, which has a red 1RAC color ([3.6/im] 

— [4.5/im] = +0.23) and is also centered on the 870 urn 
LABOCA source LESS J033302.9-274432 (SWo M = 6.7±1.3 
mjy). Since this object is undetected in the radio and MIPS 
24 ^m, it is also likely a high redshift object (zmr > 4.1). 
IBiggs et alj l|201lf ) did not find any robust or tentative coun- 
terpart for the LABOCA source. The AzTEC contours are 
significantly elongated in the east-west d irection, sug g esting 
this is a blend of more than one source. Wei ss et al.1 {2009 ) 
model this source with a fainter second component, LESS 
J033303.9-274412 (S S ro» = 5-3 ± 1.4 mjy). 

AzTEC/GS44. There is only one VLA radio source found 
within the 10.4" search radius centered on the AzTEC peak 
position. This source GS44a has a flat IRAC color, [3.6/im] 

— [4.5/im] = —0.07 and is a modest MIPS 24 [im source 
with S 2 4 M = 143.1 ± 9.3 /iJy. The 870 ^m LABOCA source 
LESS J033240.4-273802 (SWo M = 5.0 ± 1.4 mjy) is found 
within the AzTEC search radius, but its centroid is offset 
from GS44a by 12", just outside the nominal beam area 
of LABOC A. Another radio sour ce J033239. 14-273810. 5 at 
z = 0.830 (|Le Fevre et al.l 120041 ) is a better candidate for 
the LABOCA counterpart, but this radio source is located 
24" away from the AzTEC centroid, which is well outside 
the nominal search radius for the AzTEC source. We list the 
LABOCA flux in the Table, but note the large separation 
between this source and GS44a. 

AzTEC/GS45. Two red IRAC/MIPS sources (GS45b & 
GS45c) are found within 2" of the AzTEC centroid, and they 
are the most likely candidate counterparts to the AzTEC 
source, primarily by their proximity. The 870 um LABOCA 
source LESS J033218.9-273738 (SWo„ = 8.1 ± 1.2 mjy) is 
found within the AzTEC search radius, half way between 
GS45b and GS45a, the latter of which is a bright MIPS 24 
H-m source and a VLA radio continuum source with a rather 
blue IRAC color ([3.6/im] - [4.5/im] = -0.32). A higher 
resolution K-band image shows that the brighter radio peak 
is associated with a fainter component, suggesting this is 
a blended source. GS45b is 5.8" away from the LABOCA 
centroid, and it is a r obust counterpart for both AzTEC 
and LABOCA sources (|Biggs et al.ll201ll ). 

AzTEC/GS46. No radio source is found within the 13.0" 
search radius. The IRAC/MIPS source GS46a located 6.2" 
away from the AzTEC centroid is the tentative coun- 
terpart based on the combination of its red IRAC color 
([3.6/im] - [4.5/im] = +0.17; P coior = 0.255) and MIPS 
24 um emission (P 2 4 M = 0.226). The LABOCA source 
LESS J033157.2-275633 (S 87 o M = 4.8+1.4 mjy) is located 
23"from GS46a and thus also far from the AzTEC centroid. 
We list the LABOCA flux in the table but note this large 
separation. 

AzTEC/GS47. The faint VLA 1.4 GHz radio source 
GS47a (Si.4 = 43.2 ± 7.0 /iJy, P 1A = 0.105) is found 
7.6" away from the AzTEC centroid, and it is also a red 
IRAC source ([3.6/im] — [4.5/im] = +0.46) as well as a 



MIPS 24 urn source. The 870 um LABOCA source LESS 
J033208. 1-275818 (SWom = 7.3+1.2 mjy) coincides in po- 
sition with AzTEC/GS47 within 3.4" of each other, makin g 
it a robust LABOCA counterpart as well (|Biggs et al . 2011). 
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